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CHAPTER I 


ITJ TRODIT CTIOl! 

1.1 IpTRpjJUCTIOH 

In recent ye ars considerable attention lias been paid 
to the specially structured linear integer programming problems 
defined on graphs. "Good" algorithms as defined by Edmonds 
HIS]] have been reported for some of these problems. In this 
thesis we will concentrate on a generalised version of one such 
problem. In this chapter we will give a brief introduction of 
the problems of this class and later, a chapterwise brief 
summary of the remaining chapters of this thesis. 

Let G = (L,B) denote a finite, undirected and loopless 
graph where N and Ei are node and edge sets of G respectively, 
and |Hj = n, jB| = m, C is an n x 1 cost vector associated with 
the node set M of G. The problem of covering of edges by nodes, 
is to select a subset T of H, which covers all the edges of G 
and which also minimizes the cost, over all such covers of E. 

A node is said to cover all the edges of G, which are incident 
upon it. Above problem is referred to as "Weighted Edge 
Covering Problem (Cp) n . For nomenclature of this problem and 
other problems of this class we will follow Balinslci [ 7 ]] . 

CP is formulated as the following integer program by representing 
a cover T by a binary vector U, such that u = 1 if and only if 

Jr 

p 8 T : 
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CP : 


where , 


m in imi ze CU 
such that AU > e 

11 ^ = 0,1 V p G N 


A is an edge-node incidence matrix of G and e is a 


s urrimat ion ve ct or . 

Let CLP denote the linear relaxation (i.e. the relaxation of 
integrality const ra into) of Op. Dual of CLP is 


M LP : 

maximize e Y 
such that Y A < C 
Y > 0 

which is a linear relaxation of the following C-matching 
problem (MP) : 

KP : 

maximize e Y 
such that Y A < C 

Yp tl = 0,1 ¥ (p,c.l) e K. 

Another problem which is equivalent to OP, and which is also 
of great interest is the "Vertex Packing Problem (VP) ni . Vertex 
packing problem is to find a subset I of I with the highest cost 
such that no edge of G is incident upon two nodes of I. It is 
formulated as the following 0*1 programming problem : 
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VP : 

n 

maximise £ c u 
P=1 P P 

such that u + u < 1 ? (p,q) e E 

P 

u = 0 5 1 9 p S K 
P 

VP is equivalent to CP In the sense that solution of one can "be 
obtained from the other by a simple linear transformation 
U = e_U. Let VIP denote the linear relaxation of VP. 

Dual of VIP is 


minimi ze 

i 2C 




(P 

,q)eh P 




such that 


> a 

¥ p 

e n 


qSN Pa 

- P 




X pl 

> 0 

¥ (p 

CO 

✓ — \ 

O* 


Por the special case when c =1 ¥ p e H above problem together 

P 

with the integrality conditions on x * s is known as the 

P 

"Covering of Nodes by Edges" [“" V 1 • Edmonds and Johnson | 1 6 ~1 
proposed the f ollow ing general problem : 

minimise £ c sc 
(p,q}SE 

such that b 1 < AX < b g 

x p q = 0,1 ¥ (p,q) C B 

which obtains a maximum matching or minimum covering as a 
special case. 

Y/liite 58]]] discussed the problem of maximum matching and 
minimum weighted covering of fixed cardinality. 
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1 .2 SURVEY. OF EE! LAI S B PR OBLEMS 

Oonsiderable work has been done on those closely related 
problems. Berge [7 8 I] has proved the connection between minimum 
cardinality cover (when c = 1, f p e N in Op) and a maximum 
cardinality matching (when c = 1, in MP ? p c F). lie introduced 
the notion of alternating and augmenting paths and proved that 
a matching is of maximum cardinality if and only if no augmenting 
path exists. Similarly for minimum covering problems Norman and 
Babin [(45 [j defined alternating and reducing paths. Following 
are some of their main results : 


The orem _1 2 » 1 : Every irreducible cover is a minimum cover, 


'ljioorem__l_.2_._2 : If i,i^ and 1,1 g are two covers then Mg can be 
obtained by applying; a finite number of linear transformations 


on M 


1’hese notions are generalized for weighted covering and 
matching problems by Wit zg all and Zahn ("59 (] and Balinski £ 7 (] . 
Edmonds [( 14 3 has developed algorithms to search these 
augmenting and reducible paths, hitzgall and Zahn [(59]] and 
Balinski [( 7 3 presented simple alternatives to this search 
procedure, but these amplifications appear to be effective only 
for tiie problems when 0 is a constant multiple of the summation 
ve ctor . 


Vertex packing problem has recently received considerable 
attention. Properties of the convex hull of feasible solutions 
of VP have been investigated by Chavatal £12 3 » Padberg [(47 
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43 J and Hr otter 
aspects of "VP have also been studied by Balinski j” 7 3 , 

Trotter j 55 3 ? .Bellas and Samuels on jj5j . 

It has also been shown that the vertex packing problem 
and lienee the weighted edge covering problem are ” Complete" in. 
the sense of Karp f[3o I] • 

Memhauser and Trotter [_' 43 introduced the new idea of 
solving VP by a decomposition method which is not possible for 
arbitrary integer programming problems. They have also developed 
some useful optimality conditions which can be summarized as 
follows : 

The PS.SZHL ArFL’J?. C 43 Zl : 10 ^ v be a vertex packing of G. Let 
V = Li/V, then V is an optimal packing in G if and only if every 
maximal packing I C V, V(l) is an optimum pa dicing in the bipartite 
subgraph G of G, induced by I <J 7(1), where V(I) = V f| N ( I ) and 
N(l) = {p G i'i'/l ; (p,q) g L. for some q e 1 } . 

The above optimality condition then gives the following 
’local 1 sufficient condition for optimality : 

ffi&QEem 1.2 .4 |~* 43]] i If Hs an optimal packing in G, 
induced by v\Jn(Y) } then Y CV*, where V* is an optimal packing 
in G. 

By using the special structure of its constraint set, 


[^55^j . Algorithmic 


llemhauser and Trotter £ 


they have proved the following result which is useful for 
decomposing the problem into two parts. 
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1 

Theorem 1.2.5 [" 43 : Suppose X* is an optical (0, g , 1) 

valued solution of VIP and V = {p | x* = 1 1 . There exists an 
optimum packing in G that contains V. 

fho algorithmic effectiveness of the above result is 
enhanced by the fact that VliP is equivalent to VP on a bipartite 
graph G 1 and therefore can be solved by a good algorithm. This 
is a result by Kdmonds and Pulley blank, de teals of which is 
available in H'43~| . T oil owing is the restatement of this 
result which represents a natural correspondence between feasible 
solutions to VP on G 1 and feasible (o, g- , 1) valued solutions 
to VIP on G, that preserves objective values appropriately. 

T i to o r e m _ Jl . 2 _. 6 j^' 43.*] : (i) If (U, U 1 ) is a feasible solution 
to VP on G r , then U = g (U+U') is a feasible (0, g , 1) valued 
solution to VIP on G. 

(ii) If X be a feasible (0, g ,1) valued solution to VIP on G, 
then (II, U 1 ) is a feasible solution to VP on G 1 , where 

— 1 

( 1 if u = or 1 

> P * 

u = ; 

P \ 

I 0 else 

f 1 if u =1 

, ! P 

an d u = < 

P 

0 o ise . 

Furthermore U is optimal to TIP an G if and only if 
(U, U l ) is optimal to VP on G 1 . 
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Since the usefulness of the result of theorem (1.2.6) 
depends upon the number of variables having integer values 
in an optimal solution to VIP, they have also presented an 
algorithm to determine an optimal solution to VIP in which a 
maximal collection of variables are integer valued. 

Finally using the above results Kemhauser and Trotter 
[_ 43 Jj have presented an algorithm for VP. An outline of this 
algorithm is : Given a vortex packing P, an exhaustive search 
is made on V = F/V for an augmenting subset I. If one is found, 
V is redefined as (V\Jl)/V(I), and the procedure is repeated, 
il not, V is an optimal packing in G. To search an augmenting 
subset they have used the straight forward implicit enumeration 
scheme on V. Computational results given by them show that 
this algorithm performs host on the graphs of low and high 
density. Graphs of medium density are of greater difficulty 
for the procedure, not only in terms of computational time, but 
also in terms of augmentations required to attain optimality, 
Picard and ^ueyranne j_4S _] have proved that there is a unique 
maximal set of variables that are integer valued in an optimal 
solution to VHP. Using this, they concluded that the algorithm 
given by Homhauser and Trotter £ 43 ] for VIP also gives an 
optimal solution having maximum number of integer valued nodes. 

It looks that if C = e, for solving VP, it will be 
useful to solve VIP first. Pullcyblank £ 53 U has shown that 
it is not so. He has defined a 2-bicritical graph which has 
the property that the unique optimal solution U* to VHP defined 
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v. i , 

on such a graph is ^ ¥ p 0 IT. He has also shorn that 

almost all the graphs are 2-bic.r it leal. Thus for 0 = e, TO? 
is 'almost always’ of no use in solving "VP for random graphs* 

He did not give any such comment for the case, when C ^ e . 

Houck and Venmungati |~ 2 5^] also studied vertex packing 
problem for C.= e. Ihiey have also de ve loped a branch and bound 
typo oi algorithm for solving this. 'They have also suggested a 
procedure to construct a good initial solution for TP which in 
some cases is known to be optimal. Since it has already been 
shown that the vertex packing on a bipartite graph can be solved 
efficiently, they obtain this initial vertex packing for G, by 
solving a vortex packing problem on a bipartite graph G, which 
is generated from G by giving particular directions to its edges. 
They have proved that if the edges of G can be directed transi- 
tively, then tne optimal solution thus obtained is an optimal 
solution to TP on G. They have also shown that the group 
theoretical approach to integer ‘programming, when applied to 
vertex packing, yields a trivial group problem, regardless of 
the determinant of the associated basis. Using this, they have 
de ve 1 op o d s omc n i ce f at homing res ult s f or t he ir branch an d 
bound type of algorithm for TP and have also shown that it is 
possible to permanently fix certain variables at either zero or 
one and consequently to reduce the graph G. Computational results 
on this algorithm also show that the problems with density .25 
are harder than those with lower and higher densities* 
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Balas and Samuelsson £ 5 
significantly in this field. They restated the edge covering 
problem as a linear program (liPItfC) whose constraints are the 
facets of the convex hull iij. of GP (LPNG) is formulated as : 

minimize Z u 
pCB P 

such that Z u > |Q- j - 1 ¥ Q. G li 

p C Q i p 

1 d h u > d b lie F 

P P L o 

u >0, P e IT 

Jr 

whore K is the family of all cliques of G-, An inequality indexed 
by 3? is a facet of the convex hull, not associated with a clique. 
Dual formulation (LUNG) of (liPhG) is a relaxation of the edge 
matching problem 1£P, in the sense that a one-one correspondence 
can be established between the feasible solutions of LUNG and 
HP. They have developed a dual simplex type algorithm, 
introducing the notion of a dual node clique set, which allows 
identification of a partial cover associated with a dual feasible 
solution to IPNO . With the help of a partial cover, using a 
labelling method, primal infeasibility is gradually reduced, 
while integrality and dual feasibility are preserved. They 
have also reported that this method after minor modifications 
can be used to solve the weighted covering (of edges) problems. 
Till now very little computational experiences have boon obtained 
to compare these existing algorithms. 


^ have also contributed 
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1 .3 GEHERAIiI ZBD BD( 35 COVE RIN G PBOB IEM 

Generalized edge covering problem (P ) of this thesis 

>.■ 

is a generalization of the weighted edge covering problem (CP). 
In this problem, tne notion of an edge cover is extended to 
allow any positive integer in the right hand side vector instead 
of the usual entry of ones for all the edges of G. 

Ilathcmatical formulation of the generalized edge 
covering problem will be as follow : 

P : 

minimize 2 c u 

pen P P 

such that Up+u^ > r 9 (p,q.) 6 E 

Up > o, integer ?p e h . 
i . 4 smyrn^ or he ti es is 

ilotiva'tion behind the present work is obtained by 
bomhauser and Trotter's result for VP Besides the 

present chapter, this thesis contains five more chapters. 
Cnaptorwiso summary of these arc given below : 

In the second chapter we show that every basic feasible 
solution of the generalized weighted edge covering problem (P), 
is (Integer, Integer/ 2) valued, wo also prove that all the 
nodes which got some integer values in an (Integer, Intoger/2) 
valued optimal solution of IP (the linear relaxation of P) can 
retain the same values in an optimal solution of P. Thus we 
conclude that p can be decomposed into two problems : 
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(a) A linear programming problem IP, 

(b) A weighted edge covering problem on a reduced graph of G. 

Hence if a computationally good algorithm can be 
obtained to solve the weighted edge covering problem, it can 
effectively be utilized to solve this generalized edge covering 
problem too. 

In the third chapter, we have shown that an (Integer, 
Integer/2) valued optimal solution of IP can bo obtained by 
summing the optimal solutions of finite number of linear edge 
covering subprobloms on the subgraphs G^ of G. Bach such 
eubproblsm can easily be solved by solving a maximum flow 
problem on a bipartite graph generated by Maximum number 

of such subproblems are bounded by the highest value of the 
right hand side vector of p. The main result of this chapter, 
on which the above algorithm is based, is bhat, IP can be 
decomposed into two linear programming problems and the optimal 
solution of the IP is sum of the optimal solutions of these two 
du c omp o so d pro bio ms . 

In the fourth chapter, we have identified some more 
problems which can also bo solved by the above algorithm, after 
some minor modifications. This will include generalized vertex 
packing problem y transportation problem and some other similar 
problems. A natural modification of the above algorithm is 
reported for the caso when G is a bipartite graph. 

An algorithm is also developed hero, for the problem P, when the 
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underlying graph G is a tree. It is proved that in this case 
number of subproblems arc bounded by len , where k is a constant 

independent of r 's and n is the number of nodes in G. 

1 pi 

In the fifth chapter, computational results of the 
above algorithm for HP, using randomly generated test problems 
are reported. It is observed that the number of iterations are 
bounds d by O(n) . 

Last chapter of this thesis is not directly related to 
the weighted edge covering problem, but deals with the 
"Construction of optimal communication spanning trees' 1 , first 
discussed by Hu j^2 7~] • In this chapter we give some algorithms 
to construct optimal communication spanning trees, having some 
special structures. Finally wo suggest a minor modification to 
i lu ’ s alg or i thm £ 2 7 " ] , w h i c h imp r o ve s t ho c omp ut at i on al e f f i c ie n c y 
of t ho alg or it hm . 



OHAPjlER , II. 

IE COMPOSITE Oil OF' P [G;rJ 

2 . i MRo^qa ci 1 1 ok 

Generalized edge covering problem, defined in the previous 
c aap ter, is formulated as the following integer program : 

P ['G-;r ■'] 

minimize CD 
such that At J > r 

Tj > 0, Integer 

where G is an n -vector whose component c (p = l,,..,n) is 
a positive integer which represents the cost of assigning unit 
weight to the node p of the graph G. 

IT is an n-vector of v.rirbles, whose component u (p=l,...,n) 

P 

denotes the weight assigned to the node p, 

r is an m-vector of re uuirements whose component r 

y PI 

represents the requirement for the edge (p,q) of G, ¥ (p,q.) 0 IS. 
Following is the linear relaxation of P £ G^r ~] , 

LP [>;r Z\ 

minimize CU 
such t "licit AU > r 
Ii > 0 

In this chapter by exploiting the structure of P^Gjr"], 
we show that it is possible to decompose one of its optimal solution 
as the sum of following two vectors : 
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(a) A rounded down -vector obtained from the (Integer, Integer/2) 
valued optimal solution of UP G-jr . 

(b ) An optimal solution of the weighted edge covering problem 
on a subgraph of G, generated from the optimal solution of 

fp L g ?r ;j . 

Adding surplus variables U to Id? L, G-;r J , it can be 
written as : 


9 o 


minimize (Ju + 0U_ 
such that AU - lU^, = r 

>D 

U > 0 
' U 's 0 

A ifiSljLl ABOUT. ,;M blKlJOfUilE 


2.1(a) 


It is a well known result that in any basic feasible 

solution of the linear relaxation of the vertex packing problem, 

\ 

variables are (o, - 0 - s 1) valued. Details of the proof of this 

As/ 

result are available in Here we prove a similar result 

for IP £ (i;r J , by a different approach. 

lue prem J2D2 #< 1 : Values of all the variables in any basic feasible 
solution of 2.1(a) are either integer or integer/2. 

Prjoof : Let 3 be a basis of Id? Gjr ~] , when written in the form 

2.1(a). After rearranging the columns, we can partition D as 
follows : 


" -I 

Ai 

0 

Ap 
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where I is an identity matrix and 0 is a null matrix. 

For .3 to be nonsingular, it is necessary and sufficient that Ag 
be nonsingular* following two cases are possible for Ag ; 

(a) Ag is not an edge-node incidence matrix of a subgraph 
■ oi Gr * 

(b) A Q is an edge node incidence matrix of a subgraph G of G-. 

Case, a : If Ag is not ail edge -node incidence matrix, it implies 

that there is at least one row of Ag containing only one nonzero 
entry* In this case, by applying elementary operations, 
(multiplication by -1, addition and subtraction of rows) we can 

^ A. 

transform B to 3, such that 3 = R-^B 


-i 

Ag 

-“'"1 

A 

A o 

i ^ ... 


where is the product of elementary transformations of the type 

A 

described above and which are used for transforming 3 to R. Ag 
is a nonsingular edge-node incidence iiio.tr ix of a subgraph G„ of G. 

w Q 

Caoe b : ho transformation is needed and hence we take = I. 

Since Ag is a nonsingular, edge-node incidence matrix of G g 
each component of G 0 must have as many edges as nodes, lienee 
each component of G,, must be a cycle. 

Again by renumbering of nodes, we edi write 




P 

0 


A 1 

Q 1 


Qs 


# n. 
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where Q i is a cyclic matrix corresponding to the i^h component of 

A 

G 0 , ? i = 1, ...,t; and P is an upper triangular matrix with 
entries 0 5 -l and 1. corresponds to an acyclic portion of the 

graph. All the cycles of G will he of odd length, as Ag is non 
singular and determinant of an edge -node incidence matrix is 0 
or j 2 J according to as it corresponds to an even or odd cycle. 
Using elementary column operations of addition .and subtraction on 

A 

d, we get 


Brc = R 1 BK 


-I 

0 


0 

ID. 


IT, 


where I is an identity matrix and 


1. = 
x 


1 


1 - 11-1 . . .2 


Again using elementary column operations of addition and 
subtraction, we get 


■* -1 

0 

... » * 

. . — 

0 

W 1 

*'2 

♦ 

♦ 

• 


W t 


A — RR-^BL — 
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where , 



and A . = 

i 

0 


L 2 J 

Any basic solution of 2.1(a) is 

13 = (K a" 1 %r ) r 


S in co 


, ~1 

U . *— ■ 

x 
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— 1 

in 2.2(a), a is the only matrix with some fractional elements 

"1 "•‘•jL 

and these are vy 5 s only, 'luercfore each element of B r is either 

1 „ , 

sin integer or ^ integer. 

Corollary 2_.2__._1 ; If all the r ‘ s t:U: ' 0 nonnegative even numbers, 
every basic feasible solution of Lp[_G;r_j will be a feasible 


s olu t i on of P [/ G ; r __ J 
of IP jf_ G $ r ~ j will be 

2.3 3Qi.iHi PiOPERlIBS 
t 

Let U be an 
of LP j/ G;r ~\ . Let Of 


and. hence any extreme point optimal solution 
an op t imal s olu o i on of P [" G 5 r “] . 

OP TII& SOIL LION OP IP Gj r JC| 

(Integer, Integer/2) valued optimal solution 
be a subgraph of G such that for any edge 


(P»‘l) of G- } 


U 


C 

P 


+ 11 


a 


^pcL <==> (p ’ a) C G 


G may consist of one or more components. Let these components 

bo denoted by G ^ , G OJ ...,G^ . Following' w ill hold good for each 

0 

of these components ; 


"tf 

theorem 2.3.1 : For any ,i 1 < j < t , either all the u, 1 

whore p G G. are integers or all are non integers. 
j 

P r_oof ( : Let p be a node of G. and lot u b be non-integer. 

Since > 0, there will exist at least one edge (0,1) such tiiat 


t . 


'.Pi' 


Phis shows that u u is also a fraction. If there 

a 


t -h 

u p + Uf i 

is any other node p’ in G., there will be a path from this node 
to p in G ^ . Again using the above argument, it is easy to show 
that the node p’, as well as, all the nodes of the corresponding 

4 - 

path will have noninteger values in U . 


Similarly it can also 
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be shown that if in a component Cf., one node has integer value in 

d 

U , then all the other nodes of this component will also have 
integer values in 


let I = {i-^j 

G j whose nodes 


the index set 


H. 


be the index set for the components of 


s get integer values in u^. let J = { 
for the como orient o of G, whose nodes get 


• • » 5 jj. y 
t 2 

non integer 


values in L T 


let 

G_. = G ij { edge (p,q) of GjCi, such that p and q 
D <3 

both are in (>•}. 

J 

lot A be an edge -node incidence matrix for the subgraph G. 
IsnOsO can be written as : 


minimise GU 
such that Ab > r 
b > 0 

whore G and r are subvoctors of 0 and r corresponding to tne 
subgraph G. 

irrom the definition of G, it is obvious that U is an optimal 
solution for Ip£Gj r ‘ ] also* Decomposing liP £&jr 2 ] in terms 
of tho components of G we get : 

IP r "1 : minimise 1 Ql U, + 1 G. U. 

*-■ ~ J wot & k o T 3 3 


such that* 
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A, 


Ar 


L 


A, 


ti+tg 


1 


U 


TJ 


2 


t 1 +t 0 

JL (C 


t l +t 2 


,G >0 Vi = 1,2, 


whore for every i (i = l,,..,t^) A^ is the edge-node incidence 

matrix G- mid for every j , (j = l,...,t 9 ), A + ■ is an edge-node 
1 ^ C l + J 

in c iden ce mii.tr ix of G - . 

J 

therefore, restriction of to G^, whore i e I, gives an 
optimal solution for lpQ'G.*,r. ~] and restriction of U to G. 


where j G J is an optimal solution of UP LG.*, r.J. We will 

~b i "1 

denote these restrictions of U "by X , for i e 1 5 and Y , for 
j S J. Prom the construction of G., it is easy to see that 

J 

for every 3 G J , Y^ will be a feasible solution to Lp£G.;r.]J 
and hence optimal also. 


Following results exploit the structure of the constraint 
cot 2.1(a) with respect to being able to decompose its right hand 
side vector. 


Obse rva t ion 2 .5.1 : let Y* be an optimal linear solution of PI, 
where , 

PI minimize CY 

such that AY > b 
Y > 0 

and where b > 0 

Let Pi be decomposed into two problems P2 and P3, where 
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and 


minimize GY 
such that AY > b^ 
Y > 0 


P 3 


in such a 
can write 
to P2 and 
s olutions 


minimize CY 
such that AY > b g 
Y > 0 

way that b = t> 1 + b 9 and b 1 > 0 and b g > 0 and we 
Yp = Yg + Yg where Y g and Yg are feasible solutions 
P3 respectively. Then Y g and Yg will be optimal 
for P2 and P5 respectively. 


fhe converse is of course not true . 
easily be constructed to substantiate it. 

how ve decompose Y^ as follows : 

Y 5 = + {y e 


An example can 


* 


2.3(a) 


where e is the summation vector. r .. j (¥ ,j G J) can be 

d J 

decomposed, into the following two subproblems : 


IP 




minimize l _ c 11 

P es, p p 


such that u +u > r — 1 
P I '■ pq 


Pp 2 0 


? (p,4) e Gj 
?p e s, 

d 
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and 

Ll> L G j? e II 

minimize 2 _ c u 
P G Ig P P 

such that Up +u^ > 1 "t (p ,a) ® 

u >0 v p e a . 

p- D 

It is easy to see that and e are feasible solutions to 
UP G-.;r .~o and IP G- • 5 e ~’j respectively. Therefore by observation 

(2.3.1)Y^ and ^ e will also be optimal solutions for their 
re sp e ct ivu pr oblems • 

ft ov using the result that 1 ; e is optimal for IP G . je Jj , 

r*j D 

wo will prove the following s 


Ibiojoreiu J2.J3 ._2 s At least one optimal solution of PA^ will bo 
an optimal solution for the problem PAp also, whore 




pa 2 


minimise I „ c u 

p S G- P P 

ouch that u +u . > 1 ¥( p , q ) 0 j G J 

p 1 "" J 

u , unrestricted integer ¥ p G Gf . 
P 0 L 2 

minimize L _ c u 

pea. P P 
J 

sucu that u +u >1 ¥(p,ii) e G - 

1/ vi J 


u - 0,1 ¥ p e G . 

p D 


: go will prove the above theorem by constructing an 

optimal solution of PAg from an optimal solution of PA^. 
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Following can be stated for these problems : 

^ l 

(i) PA^ is bounded below as the vector ^ e is an optimal 

solution for its lineje: relaxation. 


( ii) If in the optimal solution u' ox PA-j_ all the variables are 0 

or 1, then U is feasible for Pi-ig and hence optimal. 

( iii)In case, u is not feasible for PA p there will exist at least 

one node say p., of 0 such that u <0* Phis is so, 
x 3 Pi 

be cause 5 suppose it is not true then at least one node say 


p 9 should have a value greater than one in it. In this 
case, there can not be any binding edge of PA^, incident 
upon this node p ;P and therefore the objective value can be 
improved by reducing the value of this node by one. This 
will then contradict the optimality of IJ for PA-, . A portion 
of the graph G- formed by the edges which are binding w.r.t. 
U and connected to the node p^, may look like the figure 
given below : 



Pig. 2.1(1) 
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Construct the sets S-j_ and as follows : 

= {p e G J 3 a path from to p, containing only the edges 

A 

binding 'v.r.t. U> 

B* = {(p 5 q) G G j p and q both are in and (p,q) is a 
binding edge > 

Let £i p - | = d ¥ p G S 1 2.3(b) 

Obviously d^ will be negative [’"positive "] if and only if tne 
corresponding component is negative [[positive ~| in U. 

Prom tne definition of S-j_, it is easy to see that all its positive 
valued nodes will have equal value say 'a' and all its negative 
valued nodes will have the same value, say 1 b 1 , Since a = -b+1, 
b < -1 and T J is an integer solution, all the positive valued 
nodes of S-p will have value greater than or equal to two in U. 
i’her of ore we can write 

d p = ^ fptSj 

where G {1,-13 and x is a positive non integer number greater 
than one. Having determined oho structure of 0, we will now 
construct a solution of PA 9 which has the same objective value 
as U. Tcis will be done by constructing a sequence of optimal 
solutions for PA^, starting with U. For d p (as defined above 
in 2.3(b)) following are the three possible cases : 

either (a) T cl <L = 0 
p 6 S 1 P P 



<zj» > 
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or (To) 1 c d < 0 

p e s 1 p p 

or (c) Z c d > 0 

pesj p t 

We •-/ill now show that cases(b) and (c) are not possible. 

A 

In case (a), we will be able to construct a vector U, which is 
also an optimal solution for PA^ and in which either the number 
of negative valued nodes has decreased or the value of at least 
one negative valued node has increased, without decreasing the 
value of anyother negative valued node. 

A 

A 

Case (a) : let U be defined as follows : 

tip = - d p if p e s 1 

u = u otherwise . 

P P 

It is clear from the definition of d p that the values of all the 

A 

negative u ' s will go up by at least one. Therefore either 

^ P 

A 

(a^) ; U has less number of negative valued nodes^ or 
(ag) : Values of some of the negative valued nodes of have 

gone up, while others have not gone down. 


We will later prove that U is also an optimal solution 
for PA-^. Hence for this case, proof of the theorem follows from 

A £ 

the construction of a sequence of such optimal solutions (U, U, 

A A A 

A A A 

(where ft is obtained from U in the same manner as U was 


constructed from U), for PA^ till a feasible solution for PAg is 
obtained. 
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Case^XJLL : 


2 c d < 0 

e s , p p 


==> X 2 c d < 0 


p e s. 


p p 


2 c d < 0 

p e s, p p 


( V x > o) 


■/e now construct a vector U* as follows : 


u* = u + d 
P P P 


if p e 


and u _ ^ 
P P 


It follows that 


otherwise 


(1) 2 _ c u* < 2 „ c u. 

p e & 3 » p p e o a p p 


(2) U* is a feasible solution for PA 1 . (It will be proved 
later ) . 

/V 

Thus U* contradicts the optimality of U for PA^. Hence 

A 

case (b) can never occur for U. 


Case X o ) • 

m w •»- jhif -m Brn mo hym i nmf 'J 


2 C n d p > 0 

p e s. p p 


=> 2 
p e S- 


Construct a vector U as follows 


u p = u p 


dp if p e S 1 


otherwise 
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Therefore 


2 _ 

p e <q 


3 U 

P P 


< 2 „ 
p G G .! 

J 


c 


P 


u . 

P 


This inequality together with the feasibility of IT for PA^ 

(which v/ill be proved later) contradicts the optimality of tj for 
PAj_. Tnerefore case (c) is never possible for U. Thus we have 

A 

shown that for any optimal solution U of PA^, if it is not 


feasible for PAg > following condition must hold : 


P 



d 


0 * 


Thus it is always possible to construct another optimal solution 
for PA-p which is less 'infeasible* for PAg. Constructing a 
series of such optimal solutions of PA-j_j an optimal solution of 
PAg with the same value of tho objective function can always 
be constructed. 

A 

2 * 4 FEA SIBILITY OP U, U* Ml TJ 

Por any edge (p,q) of G-., exactly one of the following is 

J 

true : 

(i) (p,q) C E 1 

(ii) (p,q) E^ and both the nodes p and q are in S^. 

(iii) ,(p,q) # E.^ and exactly one node from p and q is in S^. 

( iv) (p,q) $ E^ and neither p nor q is in S^. 

Case (j) : Prom the construction of the sets S-j_ and E^, we can 

say that one of the nodes p and q, say node p has a positive value 



28 


in U and the other node i,c. q has negative value in U, Therefore 

A A A *“* 

dp will he +1 and d^ will he ~1. In U and U, value of this 
node p is decreased by one, while value of the node q is 

increased by one. fnus u o + u q = 1, since (p,q) G E^ and 

u + u = 1 since (p,q) 6 E-i* In U*, value of node p is 

increased by one, while the value of node q is decreased by one. 

Hence 

u* + u* = 1 since (p,q) G Bp* 

C ase (ii) : (p,q) E^, this implies that there is a positive 

surplus on the edge (p,q) in tl. Since p and q both are in S-^, 

it is necessary that both must have positive values in U. All 

the positive valued nodes of S-j_ have an equal value and all the 

negative valued nodes of have an equal value. Surplus cn 

any such edge (p,q) will be _> 3, because the positive valued 

nodes are greater than or equal to two. Therefore u + u >1 

: PI" 

will be satisfied byU, U and U*. 

Case (ii i) : let p e Sp There is a positive surplus greater 
than or equal to one on this edge (pjq) } for the solution U* 

Now there are two cases ; 

(a) Up is negative and u^ is positive. 

A A 

(b) Up is positive and u^ is either positive or negative. 

A 

Gase ( a) : In U and U value of the p'k* 1 node is increased by one, 

while the value of node q is unchanged. Hence 

“ p + \ > 1 
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and u + u > 1 

ir Si ***** 


In U*, value of node p is xeduced by one, while the value of node 
q. is unchanged. Hence we still have 


Case_£bJ^ : 


U P + U 1 
Similar 


> 1. 

to case 


(iii) (a). 


■c 


Case_ ^iv]_ : In U, U* and U, values of the nodes p and q. are 
unchanged. Hence the constraint u^ + u^ > 1 is satisfied by all 
of them. 


2 . 5 IE PO MP O S I 'll OH OF l'H3 OPTIMAL SpiIJTM 

Using the above results, we will now prove the following 
theorem : 

theorem 2.5,1 : (V" 1 + (p) is an optimal solution of 

(¥ j G J), where Cr is an optimal weighted edge covering for 

G- and P is, as defined in section 2.5. 

J 

Pro of : let Ip be an optimal solution for P [p-pr . “) an I I e "^ 
d = - Y^ where Y^ is an optimal solution for Ip£G..jr_.[] as 

defined in (2.3). Y^ can be written as 

Y* 3 = vi + \ e - 2.5(a) 

All the half spaces (u +u > 1) (¥(p,q) e G. ) passing through 

p y. * j 

l i 

the point g e, generate a cone (with vertex g e ), which is the 

solution set for IP j[3 G- I e Z without nonnegativity restrictions 

imposed on its variables. All the half spaces u + u > r 
1 _ . P <1 — POL 

(V (p,q) e G.. ) passing through the point Y 3 generate a cone Z g 
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with vertex at Y''. Zg is the solution set for JLiP G ^ 5 r? when 

nonnegativity restrictions are ignored. Since r > 1 ¥(p,q) e G. , 

Ph j 

1 

Zg is contained in Z-^. Hence d + ^ e is contained in Z^ and 
is also integer valued. This may or may not he feasible for 
Ipj^Gve^j, but is definitely feasible for PA-^. 

By theorem 2.3.1 


Z _ (cl + 5) > 

P e ^ 


2 „ c n 0 d 
p G Sj P P 


= => Z _ c d 

p e g. p p 
D 


> 


2 - c p (Op ” 2 ) 


P 6 


2.5(b) 


Now we want to show that W 3 is no better than V 3 + q 3 . 
let 

Z _ c < Z . c (v 3 + O'b 

p 6 G-j P P P G G i P p p 
J J 


=> 2- o„(yNdp) < Z _ (t 3 + ol 


p e s, S '-B 

D 


p e g. p 

ei 


==> z _ o a 
p e Sj p p 


< 1 - (o n 3 -l, 

P 6 1 


'P 2 


which contradicts the inequality 2.5(b) 

optimal solution to P £ G .;r . J , 

J J 

It is obvious that 


IT 


( U X 1 , U (V 3 + O. 5 )) 

i e i j e J 


Hence V 3 + 0 3 is an 


is an optimal solution, for pQG$r_J. 
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It is also easy to sec that the requirements of all the 
other constrfints corresponding to the edges which are not 
present in QC-jr U aJ0e met 13 .T U*. Hence U* is an optimal 
solution for P [_ G^r ^ . 

Thus we have reduced solving tnis integer programming 
problem P L&;r ^ , into two parts. In the first part, we solve 
a linear programming problem, which is a linear relaxation of 
the given integer problem and in the second part, w e solve a 
we ighte d e dge c o ver ing p r oble m . 

Since the optimal linear solution would quite often 
reduce the graph G- into several components, above weighted edge 
covering problem can be solved independently on each subgraph. 
Hence any efficient algorithm [_ 5 , 25 , 36 , 43 for Ike 

weighted edge covering problem can be used to solve this 
generalized edge covering problem too. 



CHAP HR HI 


ALGORI'lffi FOR LP [C^rJ 

3.1 MROJJU'CIM, 

In this chapter we discuss the problem ItP [ G;rJ which 
is defined in chapter two. Although Simplex method can be used 
for solving it, we here present another exact algorithm, which 
is also iterative in nature. Inis algorithm is specially 
designed for exploiting the structure of the constraint set of 
LP [Hjr ~] , In each iteration we solve only a maximum flow 
problem on some symmetr ic bipartite graph, generated by a 
subgraph of G. 

3 .2 NOIATIOLiS 

I or the i’k* 1 iteration where i > 1 we define 

max l(i) = max {r (i-1) > where, r (i-1) is a 
( P ,q)eG pq 

component of the updated requirement vector 

r ( i-1) defined in section 3.3. 

max 2(i) = max { r ( i-1) Jr ( i-1) ^ max 1 (i)> 
(p,cOe& pcl Pli 

r ^ = max 1 ( i ) - max 2 ( i ) 

s^ = Qnax l(i) - max 2(i) [] 

where Qmax l(i) - max 2(i)^] is the smallest integer greater 
than or equal to (max 1( i) - max 2(i)) 

N G (p) = {Q.J 3 (p,q) G G> i.e., it is a set of all the 
neighbouring nodes of p in G. 
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G. 


n 

U* 



r(i) 


A subgraph 
consisting 

re quire ment 


of G generated for the i^- 1 iteration, 
of all those edges of G, for which the 
up late cl after ( i-1 ) ^ 1 iter at ion , is 


max 1 ( i ) * 


A bipartite graph generated by G^, whose node set 
is IT-, 21,. where IT-. = IT 0 and IT-, is the node set 
of G-, and (p,q) e G(G i ) if and only if (p,q) e G^ . 
Requirement vector, updated after i^ 1 iteration, 
formula for updating r(i) is given in 3.3. 


r(o) ; Given requirement vector. 


3.3 OUi U III OF SHE ALGOHIiill Ai'i'D SOLE; RElLA'lEI) HSSU113 


let U(l) be an optimal solution of U? [( G-^r^e ~J f° r 
the first iteration, he update the requirement vector r(o) with 
respect to the solution 11(1) . r(l), the requirement vector, 

updated after first iteration, is obtained as follows : 

r pci (1) =r pa (0) -U P (1) - U <1 (1) 

if r pq(°) - U p^ 1 ^ + ^C 1 ) 

= 0 otherwise. 

ye now work with IP £ G s °r (1) [J and repeat the whole procedure 
till all the requirements have been met. At any intermediate 
iteration i, we calculate r^ and find G^ from the requirement 
vector r(i-l). Solution of IP L^jy r i e IJ can chained by 
first solving the problem ia?['G^ 5 e]] a11 ^ ‘kken multiplying it 
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"by r^« Optimal solution of DP £ G.jje ~J will "be obtained by 
solving a maximum flow problem on G(G^). Proof of this is 
based on the results obtained by llemhauser and Trotter £43_]» 
By using the transformation, (described in cuaoter one) for 
reducing a vertex packing problem to an edge covering problem, 
the theorem by llemhauser and Trotter £43 [] can be adopted 
for an edge covering problem in tiie following form! 

Th eorem 5 ,5 .1 : (i) If is a feasible solution of 

P C G ( G i) ? e H J then U = is feasible (0, -g , 1) valued 

s olut ion to IP H G- • 5 e I] * 

— ^ 1 p* 

( ii) If U is a feasible ( 0 , rr ,1) valued solution to HP Q ^ ^ I e | 

then (U ,U ) is a feasible solution to . P £ G(G- ) je j where, 

1 if u . = i or 1 

1 j - ^ 

^ j 0 otherwise 

(1 if u • = 1 

2 J 

and = '■ 

otherwise 

Furthermore U is an optimal solution toIp£_’G^;e”| if 
and only if (U' 1 ',U^) is optimal to P £ G(G^) je . 

After establishing this one-one correspondence between 
solutions of and P[&(G^);e ], we obtain the optimal 

solution of F n^"(^i)5 e H by solving the dual of IP £G(G^) l e U 
which is a maximum flow problem on G(G^), 

Dual of IP £ G(G^) je [] is : 



maximize 2 „ x 

(p,l) e 8(Gp P8 

such that 2 x „ < c ¥ p G Hi 
»,t P M. P -*• 


qGI'T, 


2 x < c 

pgh/^ q - 


q e iir 


x pq _ > o ¥(p,q) e G(G i ). 


Let X be an optimal solution for DLp £ G( G^) $ e ”j . It 

is snown in Appendix A that if X is not a basic solution, even 

then, we can directly construct a (0,1) valued feasible solution 
12 - __ 

(U ,U ) of LP L G(Gj_) « e 2] which satisfies all the complementary 

slackness conditions with respect to X, and hence the solution 

thus obtained will be an optimal solution for LP []Gr(Gr. )je "”] , 

Thus a (0, g- , 1) valued optimal solution of LPj^Gjje J is 

obtained. Let U(i) denote an optimal solution of LP £Ck ;r^ e J , 

Here U(i) = r^ U, where U is solution of LPLG-^je”] and is 

1 2 

available in terms of (U ,U ). 7/e will prove that TJ* is an 

optimal solution where for some k 


U* = U( 1) + U(2) + ... + U(k) 


7/e will consider the following possible cases s 

(1) All r^*s are integers 

(2) Some r^’s are integers and some are non integers. 

7/e will first show that the proposed algorithm gives 
an optimal solution for case (1), in a finite number of iterations. 
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Outline of the_ proo f for Opti m ality a nd Fini teness : 

r * s are given to be positive integers, and at every 
pq b e 5 

iteration i, some of them are being updated by 'in 1 which is 
also a positive integer, 'therefore all the requirements would 
be met after a finite number of iterations • Let k be the total 
number of iterations needed, he consider now the following 
two problems : 

* Subproblem being solved at tne i^h iteration,- 

i.e . 

L ^ G ij r j e li 

minimize 2 c u 

p e G i p p 

such that u +u > r., ¥ (p,q) e &• 

ir “X ^ 

u p > o, t P e g. 

1 Lhe linear programming problem, with the right 
hand side vector denoting the updated requirements after i^ 11 
iteration, Inis problem will be denoted by 

^ C G $ r (i) II : 


minimize 

2 c 

u 


P 

G G- p 

p 


such that 

U P +U 1 * 

r P 4 (1) 

¥(p,q) 0 C- 


u d £ 

0, 

¥ p e g 


Let U B (i) = U( i+1) + U(i+2) + 
Since U B (i-l) = U(i) + U B (i) 
U B (o) = U*. 


• • • 


+ TJ(k) 


3.3(a) 
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B 

It can easily "be verified, that U (i) is a feasible 
solution for lip G;r( i) j. U 3 (i) +U(i) will be shown to be an 
optimal solution of IP [^G|r( i— 1) J , ¥ i = l,2,...,k. let G- 1 be 
a subgraph of G containing all the nodes of G- and all its edges 
which are binding with respect to U*. 

for a given i, where 1 < i < k, let 

S = { P J u p ( i ) > 0, but u 3 (i) = 0 } 

Y/e now prove some of the properties that are satisfied by U(i) 

B 

and U ( i) . 

(Pi) : If (p,q) e G^ and it is binding with respect to U(i) i.e., 
U p (i) + u q (i) = r i 

then the updated requirement on this edge will be the 
highest. In other words, 

r pq( i) = max 1 (i + 1). ' 

Proof : Por any edge (p ! ,q‘) following three cases are possible : 

(a) (p ! ,<!’ ) $ G jL 

(b) (p'jl 1 ) S Gq and it is binding w.r.t. U(i) 

( c ) (p'jl 1 ) G Gq and it is a non-binding edge w.r.t, U(i). 

Case.. ( a) : (p ' » 1 T ) ^ Gg- =:= > r !n ,(i-l) < max 2(i). Since 

x P 4. 

Up t ( i) and u ,(i) > 0 

r ptq«(i) < max 2(i). 

CaseJ^bJ : (p',q’) e G j[ ==> r pI ^,(i~l) = max l(i) and since 

(p r ,q ! ) is binding w.r.t. U(i), 
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**• r p iq_!(i) = max {r pIqt (i-l) ~ u pI (i) ~ u p! (i), 0 > 

= max {max 1(1) - max l(i) + max 2(1), o> 

= max 2(1) • • max 2(i) > 0. 

Gas.® 1?J. • (PSI 1 ) 6 and let u p , ( 1) + , ( 1) = r + e 

for some e > 0 

Since r pI ^ f (i) = max { r p , ql ( 1-1) -u p ,(i) -u pl (i), o> 

r p >q' ( 1) = ma:x {max 1 ( i ) - ^i - G, 0 > 

= max { max 2 ( i ) - G , 0 > 

< max 2(1). 

Hence we conclude that 

max l(i+l) = max 2(i). 

( 12) : lor every iteration 1, all tne edges of G corresponding 

to the positive components of r(i) and which are also 

T> 

binding with respect to If (i), will also be present in 

G 1 . 

Proof : Let (p,q) be such an edge. By definition 

r pd (1) = l?q " 1 ± (U p (3) + V 3)) 

Therefore if (p,q) is a binding edge of U B (i) and r (i) > o, 

P Q. 

then 

+ + ^ (yj) + ys)) = y- 

But by equation 3.3(a), left hand side of this equation is 
u* + u*. Hence (p,q) G G 1 . 
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(P3J_ : If a node p is in the set S, then all the edges of 

Ip["G^jr^ e J which are incident upon p and which are 
also binding with respect to U(i) will be present in 

G 1 . 

Proofs let (p } pi) , (pjpg); .. m(PjP s ) be a11 bhe binding edges of 
LP^G^r^ e , incident upon p e S. 

In proving (Pi), we have shown that 

max l(i+l) = max 2(i) for any i such that 1 < i < k. 
‘therefore we can write 


max l(i+l) = ^j_ +1 + r i+2 + ••• + r j c 3.3(b) 

3 

Since U (i) is feasible for IP G;r( i) ]]] it will satisfy the 
f oil ow ing c onstr aints 


u p(i) + u p .C 1 ) 2 r pp .(i) * 3 = l,... 5 s 

3 3 

But by property pi, r (i) = max 1( i+1) ¥ j = 1,2,'. . ,,s 

pp 3 

B n 

and since p e S, u^(i) = 0, therefore Up_.(i) > max l(i+l) 
¥ j = We will now show that u-^ ( i) = max l(i+l) 


f 3 = l,2,...,s. Lot (I) > max l(i+l) for some j o ,l < < s, 

^o 

then by 3.3(b), there exists at least one t, ( i+1 < t < k), such 
that 


u p . (t) > r t 3.3(c) 

3 o 

Since c > 0, it will be possible to reduce the value of CU(t) 

P D 0 

But this will contradict the optimality of U(t) for IP [[G^jr^ e J . 
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Hence uT 3 (i) = max l(i+l) ¥ j = l,2,...,s. How 
P 3 


l* -i- u* 
P P 


x 
2 

z-l p 


1 

2 

2 = 1 P 


u-' -!- U* F ( 2_ U n (z) + 2 u n (z)) + (Up(t) + (i)) 


ID 

X 

P 


— (x •— iicizc 1( i+1) ) x- max l(i+l) 

pp d 


= r 


‘PP 


Hence (p ,p -) G G 1 ¥ 3 = 1,2, 

£P4J^ : Let be a node adjacent to a node p e S, such that 

(Pl?p) Is an edge of G-^, and binding v-.ith respect to 
U(i), then no node adjacent to p^ in G 1 can take a 
positive value in TJ^(i). 

P£oof ; Let ( p -j_ ? q.^ ) 8 G* . Since (p,p-j_) is a binding edge 
with respect to U(i), by property pi. 

^pp (i) = max 1( i+1) 

As p 8 S, by property (P3), 

B 

u (i) = max l(i+l). 

P 1 

B B 

We will no w snow that u (i) = 0* Letu (i) > 0 

ll Q-i 

u® (i) + u® (i) > max l(i+l) 3.3(d) 

P 1 ei l 

I 

Since r p q - 2 (u (j) + u ( j)) < max l(i+l). 

L j=l p l 4l 

Therefore 

r Piqi - + «4 X ) < «*+*> ~ O^U) + 
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But by 3.3(d), right hand side of the above inequality is a 
no g t r ve numbe r , 


- p l c il 


< ud 


Pl 


+ u 


fhis contradicts the assumption tnat (p^jlq) Q G* . 

.IQ nee uP (i) = 0» 
q l 

(Pj5j( °° Bet (Igj do > • • > 136 a subset of positive valued nodes 

in U(i). If v/e get a feasible solution of BP G^jr ^ e ~p \ , 

from U(i), by reducing the values of each of the nodes 

1 1 

q.1 , q.g , . . . , q^ by g- and increasing by g , only the values 
of noues in B(q^, q^, . . . ,q.j_) , where 

B(ci 1? q 2 , ...,q t ) = {p|(p»< 3.j) is a Binding' edge w.r.t 

U (. i ) f or s ome 3 = 1,2 , ...,t> 

then c +»».+c < I c. 

q l q t " P S B (q 1 , p 

Pr o of j Follows immediately from the optimality of U(i) for 
LP e ~] . 

The ore m_ 5 .5 »_5 : U(i) + IJ B ( i) is an optimal solution to 

LP £ G’ jr 1 ( i-1) P\ if U B (i) is an optimal solution to LP £ G’ jr * ( i) “j , 
where r'(i) is a subvector of r(i) corresponding to the edge 
set of G 1 * 


Proof : It will be proved by contradiction. Let T be an 

optimal solution of 13? [_*G‘ 1 ( i-1) such that 

CY > C(U( i) + U B (i)). 
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Since c p >0; ? p G P'j there will exist at least one node p 0 , 
for which 


u (i) + u.- 8 (i) > 


Since Y is an optimal solution of a linear programming problem, 
we can select Y to be a basic solution, therefore by theorem (2.2.1), 
Y will have integer, integer/2 values for its components. 

Ye partition the node set of G 1 , into the following subsets ; 

Si = (P | p 8 S and y p < u^(i) + u p ( i) > 

= {p | p is connected to a node of S^ by a binding edge 
of e I! mid y p > u.p(i) + u^(i)> 

s 2 = { P I u p( i) > 0 and y p < u^/i) + u p (i)> 

Sg = {p j p is connected to a node of Sg by a binding edge of 

problem and y > u (i) + u^(i)} 

Si = {p ) P / sj U Sg U S:f U Sg, and y p > u p (i) + u p (i) } 

Sg =' {P 1 y p = U p (i) + 1^(1) } 

Thus ¥ p e sj y Sg, 

y p < u Q (i) + u-p ( i ) , 

SjXI d. t * I T \ 1 

¥ p e sj \j Sg \J s 1 

y P > Y (i) + 

vie now prove some properties of these sets : 
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(a) 

(b) 


(o) 


(d) 

(i) 

and( ii) 

Case (i) 
edge (p,< 


S 1 f\ S 2 ~ * 

S i (\ S j = * ¥ *•» 3 e U,2 > 

These follow by the definitions of S^, S^, Sg and Sg . 

s i n = <■ 

Proof is by contradiction. 

let s[ H S 2 ^ ^ 5 fmd let P 0 e S 1 H s 2 

p Q G S^ ==> p is connected to a node of S, by a binding 
edge of iP^G^r^ e . 

P 0 6 Sg ==> p Q should be adjacent to a node q Q in G* 

such that ui 5 (i) > 0. But by property (P4), no 
cl o 

adjacent node of p Q in G T can be positive valued in 
U S (i), therefore Sg f'l Sg = . 

It is easy to see that every constraint of IP £ G ’ ,r 1 ( i-1) 
is binding with respect to U(i) + TJ^(i). Therefore, 
there will be no edge of G' which connects two nodes 
of S* Vj Sg yj S± and is also binding with respect to X. 


Similarly there can not be any edge (p-j_»q-j_) of 
such that p-j_ e Si \J Sg \J S 1 and q^ 0 Sg. 

We will now show that 

~ ^D c n + ^I c a — ® 

p e p q e s£ q 


G 1 


^ 1 °D + ^ I c a 0 * 

p e s£ p q e sj ^ 


: p 6 sj ==> p 0 S, therefore by property (P3), every 
.) which is binding with respect to U(i), will be 
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present in G 1 and therefore it will be binding in iiP [" G,r ( i) 
also. Since 

5p < V 1 ) + u p (i) 

an d v + y > r ( i-1 ) 

P 1 — pqp ' 

= u p (i) + u q (i) + Up(i) -i- u^(i) 
therefore > u^(i) + u®( i) * 

Hence q e S^. 

Thus we can say that for ever y node p e S^, all the 

neighbouring nodes q. of G-, such that (p,q.) is binding with 

J J 3 

respect to U(i), will be present in S^. 

Hence by property (P5) 


2 -p. c -f l j c 0 

p G Sf P 9 8S 1 ! 1_ 

S.a#®JLi.U. ! Pi e s 2 ==> u p t 1 ) > °- 


By property (P2), all the edges, incident upon and 

TO 

binding with respect to U (i) are present in G } . let (p 1 ,q 1 ) 

be one such edge. As in case (i) v/e can show that q 1 e 3g. 

B 

Since U (i) is assumed to be an optimal solution of IP j^G 1 jr 1 ( i) 3 


^ D c n + ^ I c n 2 0 

P e S£ P q e si a 


Hence, we can say that 


1 -p, ; . - n c - I j , • t 0 < 0 

p e s? \j s5 p q e sf Ij s* ^ “ 
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Let us now construct a vector Y, from Y Toy making the following 
additions and subtractions % 


(1) Decrease the y 1 s by a constant d for all the nodes p 

P 

of S:f\J Sg where d < 

(2) Increase the y^ 1 s by d for all the nodes p of U Sg. 

By property (P3), it is easy to see that Y is a 

feasible solution of 13? £ G ! jr( i—1) ^ and 


~CY + OY = d 


p e s? t J s 5 ° p q e s? \ J si q s e s 


c - 2 


1 


1 


Since the right hand side of the above equality is a negative 

number, a contradiction to the optimality of Y, is obtained. 

J3 

Hence U(i) + u (i) is an optimal solution to LP L G-‘ jr 1 ( i-1) ^] , 
B 

provided U (i) is an optimal solution for IP £g* jr ! ( i) • 

[Theore m 5.5.4 j U* = U(l) + U(2) + ... + U(k) is an optimal 
solution for LP^G-jr^* 



Proof : y/e will first show that U* is an optimal solution to 
I® C&'iX* □, where r ! is a subvector of r, corresponding to the 
edges of G* . We, will prove this theorem recursively. 

St ep 1 : U^(k-l) is an optimal solution for 13? £ G 1 ;r 1 (k~l) ^ : 

Since 1c is the last iteration, all the left over 
requirements would be met by U(k). 

Problem solved at the k^* 1 iteration is 
minimize 2 o u 

p e djj. p p 

such that u p +u q > max V (p,q) e G-^ 


•^p e G k 
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Since all the edges of this problem, which are binding 
with respect to U(k), are present in S', U(k) is also optimal 
for liP ^G-’ jr 1 (k-1) 2] . hence U JJ (k— 1) = U(k) is optimal for 

IiP C G* ;r »(*-!) 3 • 

Step_2 o By theorem (3.3.3), U JJ (lr-2) = U(k-l) + U B (k-l) is 
optimal for BP [2 G- 1 ;r 1 (k~2) 2] j because U B (k— 1) is optimal for 
BP £G« jr * ( k— 1 ) 23 . 

S tep 5 ; Let U B (k-i) be an optimal solution for BP Q’ G 1 ,*r 1 (k-i) . 

By theorem (3.3.3) U B (k-i~l) = TJ(k-i) + U B (k-i) is optimal for 
hP Qg r ?r 1 (k-i-1) 2] • Hence we conclude that U B (o) = U* is an 
optimal solution for BP[2G I ;r ! 2]* 

From the construction of Bp[]G$r(i)21 for i = 1,2, ...,1c, 
it is easy to see that U* is a feasible solution to BP [^Gjr 21 
also. Since BP^G^r 1 H has fewer constraints,!!* is also an 
optimal solution for hP [2&;r 2] • 

□ 

3.4 MPIHIPBnBSS OP THS PROPOSSID ALGORITHM, WHEN SOUS r.«s 
' “* ~ — 

It is easy to see that for ary iteration i of the above 
proposed algorithm where (2 < i < lo), number of nodes in G_ j . is 
greater than or equal to the number of nodes in G^ but it is 
quite possible that the number of iterations with the same set 
of nodes may not be finite if some r^s become fractions. Por 
example consider the BP problem on the following graph : 

0 
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2 


2 



1'ig . 3.4(1) 


(Number beside an edge represents its requirement) 
Let c =l ¥p=lj2 ? ...,6. 

±r 

Iter ation 1 : max 1(1) =3 

max 2(1) = 2 
r ^ =1 

Therefore G-^ is 


c* 



An optimal solution for 13? LG-^j© H is 
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An optimal solution for IP £ ej is 



max 1(3) = 3/2 

max 2(3) = 5/4 
r = 1 
^ 3 4- 

Optimal solution for IP e £] is 


i 


l 
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U(3) 


0 


! 1 



Thus for any large number i Q , optimal solution of i Q ik iter ait ion 
w ill be 

| o 1 

; i I 


UUJ = 


and r ( i ) f 0 . 


i 1 

! ?< 

0 

1 


i. 


2 

0 


Thus if some r^s are not 
not terminate in a finite number 


integers, the algorithm may 
of iterations. 


Following modification of algorithm (3.3), will ensure 
its termination to an optimal solution of IiPL^D in a finite 
number of iterations. 


Algorithm (3 .4) i 

Here we use Sj, in place of r ^ in the algorithm 3.3, where 
has been defined in section (3 .2). All the other steps, will 
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be exactly the same, as in algorithm (3,3). 

P toit eness^ of Al gorith m j[3j,4J_ : 

Initially all .r ’ s are given to be integer numbers. 

Since at every iteration 'i T some of the r l s getting 
reduced by an integer number Sp, all the requirements would be 
met in a finite number of iterations. Hence this algorithm 
(3.4) will terminate after finite number of iterations. 

let here also U(i) denote an optimal solution of the 
subproblem IP j/G-pi s^ejof i^ 1 iteration. By definition of the 
subproblems for i = 1,2, it is easy to see that U* is a 

feasible solution to HP j^Gjr /j . 

We will now show that property Pi (proved for algorithm 
(3.3)) will hold in this modified algorithm (3.4) also. 

Proof % Since r *s are given to be integers, r. can become a 

— p Cj_ ^ 7 1 

fraction only if solution of some preceding subproblem was non- 
integer. let i Q denote the first iteration for which U(I 0 ) is 
(integer, integer/2) valued. Thus ¥ j < i 0 , max l(j) and max 2(j) 

both are integers and s . = r ., therefore Pi holds upto i 

J J ^ 

iteration. Hence max l(i 0 +l) is integer and max l(i 0 +l) = 
max 2(i Q ). if max 2(i Q +l) is also integer, Sp +1 will again 

0 , i 

he same as r. , and therefore Pi will hold at ( i +l) tn 

O 

Iteration also. We now consider the case when max 2(i Q +l) is 

not integer. Here in the modified algorithm (3.4), for (i 0 +l)"kk 

iteration. We replace r- , 1 by s. , 1 , 

x o x o 

where Sp +1 = max l(i +1) - max 2(i +1) + ^ 
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We will now show that Pi holds for (i +1)^ iteration also; i.e. 


max l(i 0 +2) = max 2(i o +l) - ^ , 

and max 2(i Q +l) - ^ is the requirement on all the binding edges 

of G. In r(i +1), updated renuirements for the binding edges 

x o ° 


of Gr • -1 8X6 

1 o +1 


= max l(i +1) 


Ji O +1 


1 


= max 2 ( i Q +l ) - g 


It is easy to see that in r(i +1), there can not be any 
component having value greater than max 2(i +1). [Therefore it 
will be sufficient to sho w that 

Lemma (3.4.1) : max 2(i.+l) ~ i = max {r (i +1)>. 

8 ( P)4 )e g 0 . 

Proof : Edges of LP 0 ;s- e/] which are binding w.r.t. U(i ), 

1 o 1 o 

either form one connected subgraph of G. or two or more 

i° 

components. As proved in theorem (2.3.1) either all the nodes 

of a component are integer valued in U(i ) or all are integer/2 

valued. Let <J(i ) be the index set of all the above components 

having noninteger valued nodes and let F^ ; j 6 J(i Q ), denote 

these fractional components. From the definition of i , it is 

clear that every node p of 0 F . is fractional valued 

3 6 J(i 0 ) 3 

in U(l) + U(2) + ... + U(i 0 ). Hence all the edges of G, having 

fractional requirements in r(i ) must be incident upon a node 

» i , , 

of \J F.. Since pi holds for i " cn iteration, G. 

3 e J(i 0 ) 3 0 V 1 


will 
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contain all the edges of P . and all those other edges 

3 e J(i 0 ) 3 

(p,q) of G for which rp^(i Q ) = max l(i Q +l). Vie will now prove 
the following lemma ; 

Lem nia (3.4.2) ; Ho new edge of G (i.e. edge which is not in 
G. ) incident upon a node of \J P . can he present in 


3 e J(i 0 ) 


% +1 ‘ 


Pr oof : let (p^jPg) 821 e< iS e which was not in , and is 

i t 

incident upon a node of P-. 

3 e J(i 0 ) D 

Since (p 1 ,p 2 ) $ & i 


r p!P 2 (i o“^ ^' max g CV- 

Srnce r pip ^(i 0 ) = max {r Pl p 2 ( i 0 -l)-u p ^(i 0 )^u p ^( i Q ) , 0 > at least 


one node from p-, and p Q will he in \J P therefore 

1 3 6 J(i Q ) 3 

u ( i ) + u ( i ) > 

n . v Cs ' \ — r 


Pi 


Pg ° ; “ 2 * 


Hence 


r Pl P z (t <P i max s < x 0 ) “5 

< max l(i +1) 


• (Pq>P 2 ) $ -j-p 1 


We now continue with the proof of lemma (3..4*«1). As a 

result of lemma (3.4.2) we conclude that all the fractional 

valued components j 0 J(i ) of G. are also present in 

J ^ 
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G- as its components. Since for solving the subproblem 
x oJ 

XiP l_.Gr. ; e J 5 we first solve H+’G. je ~] and then multiply 
1 o 0 1 o 

the solution vector by s. , restriction of ■ U ( i +1) to 

1 _ G ■ . -i 0 

0 1+1 

0 

P- is same as the restriction of -~-u(i ) to P . . Therefore 
3 s . v o' i 

V p e 


3 e J(i 0 ) 


p . 


s - 


u (i +1) 
P^ o 1 


V 1 


3 . 4 ( 1 ) 


Hence the requirement on any edge incident upon a node of 

vJ Id gets updated by atleast s- -,/2 in r(i +1). 

3 S J(i Q ) J V 0 

How let (p,q) be an edge of G 


(a) (PjQ.) 6 Gr. 


i +1 
o 


==> r pp (i o) = max 1 ( i n +1 ) 


0 


therefore r pq (i 0 +l) = max (max l(i Q +l) - u D (i Q +l) - u n (i 0 +l),o > 


P 0 


q' o 


< max {max l(i +1) - s. 0 } 

1 o 

1 

< max 2( i Q +l) - g- . 


(b) Let (p,q) G CJ- - G. - . 

o 


Since (p,q) # G. 


i +1 
o 


==> r p i o^ i max 2 (VI) 

Again two cases are possible 


r p4 (i 0 ) < max B(i 0 +1) 


either (i) 
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or (ii) r (i ) = max 2(1+1) 
v ' pq v o' '■o' 

Case (i) r (i ) < max 2(1+1) 
v ' pq v o' ^ o ' 

==> r (i +1) < max 2(1+1) - i . 
p q v o ' ~ v o ' 2 

Case (11) r pc[^o^ = ma:c 2(i Q +l) = a P rac "ki° na i value, 
lien by 3.4(1) 

^1+1 

r ( 1+1) b max 2(i +1) — 
pq^ o ' — v o ' 2 

1 

< max 2( i Q +l) - r> 

Hence 

max 2(1+1) - i = max {r (i +1)> 

0 2 (P,!) e G 0 

Thus max l(i Q +2) = max 2(i 0 .+l) - ^ and Pi is true for (i^l)^* 1 
iteration , 

In short we can say that pi holds upto (i +l)kk iteration 
because of the following property ; 

: Por iteration i 0 all the edges having fractional 

requirements in r(i ) are incident upon a node of V/ P - 

d e J(i 0 ) 3 

1 

and therefore their requirements get updated by at least ^ in 
U(i 0 +1). 

Now to prove that Pi holds for other iterations j > i 0 +l also, 
we will sho w that every requirement vector r(j), j > i Q +l has 
the property o. Consider the (i +2)"^ iteration. Two cases are 
p oss ible : 
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(a) max 2(i +2) is integer. In this case, 

s^ + 2 - r j_ + 2 = nia2: l(i Q +2) - max 2(i 0 +2) 

therefore pi holds for ( i +2)"^ iteration. 

v o ' 

(h) max 2( i Q +2) is non integer, let (bqjPg) be ^ edge of G 
having requirement equal to max 2(i 0 +2) in r(i Q +l). 
Consider the following cases : 


(i) s i +1 is e7en 

o 

(ii) .q is odd. 
x o + 

Ca se ( i) : In this case U(i +1) will be integer valued. 

Therefore only the nodes of 0 P . will have non integer 

j e J(i 0 ) 3 

values in TJ(1) + U(2) + ... + TJ(i +1). Hence any edge of G which 
has noninteger requirement in r(i Q +l), is incident upon a node 
of 0 I’ ., where J(i +1) = J(i ). 

j e J(i 0 +i) 3 

Hence property o holds for ( i^+l)"^ iteration. By using lemma 3.4.2 


it can be proved that all the components IV, j e J(i Q +l) will be 
present in G^ + g as its components. Therefore in U(i +2), every 


node of 


( j i.+2 1 

v.' P . gets a value — — — - > 

3 e J( i 0 +l) 3 8-8 


and hence 


r PlP2 (i ° +1) ^ maX 2(i o +1 ) ~ 2 * 


Thus Pi holds for (i +2)”kk iteration also. With the help of 
arguments made for ( i +1)^ iteration, we can show that property 
0 will hold for all the iterations i > (i +1) also. Hence 
property pi will continue to hold. 
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£.P, s j Liiy-Jt : s j_ + i is odd, therefore — is fraction. All 

\ \ i +1 

tae nodes of p will get the value in U(i +1)* 

ieJ(i 0 ) J 2 0 

Therefore all the nodes of P . vvill have integer values 

3 e J(i 0 ) 3 

in U(l) + ... +U(i 0 +1), and J(i Q +l) f\ J(i Q ) will he empty. 

Hence it can he stated that all the edges having fractional 

requirements In r(i +1) are incident unon a node of \J F.. 

3 e J(i 0 ) 3 

Thus property o again holds for (i o +l) t]a iteration and therefore 
PI also holds for (i Q +2)"^^ iteration. We can now treat (i Q +3) r< i 
iteration as the (i Q +l)st iteration and repeat the arguments to 
show that PI holds, for every i (i = 1,2, . . . ,k) .. 

□ 

All the other four properties P2, P3, P4 and P5 follow 
from pl, therefore thejf hold in this case also. Hence theorem 
(3.3.3) is also true for the case when some r^’s are fractions. 

The orem JJ> .5 .4) : 

U* = U(l) + ... + U(k) is an optimal solution 

f or IP C G ; r "2 . 

Pyool i Feasibility and optimality of U* can he proved on 
the same lines as in the proof for IT* in algorithm (3.3). 
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3'. 5 

Algorithm (3.4) to solve Lp [GqrJ 
Input G £. (I'-iyg) 

requirements {r > ? ? (p,g) 6 E 

P A 

costs Cp 5 T f p e E’. 

Step_ P 1 Initialize index count : i ^ 0 

r ( i) r 
P c i ; <- p q 

B(i) - El 

Step^jL : i <- i+1 

max l(i) - max {r (i~l)} 

(p,i) e s p c1 

B(i) - {(Pjl) | r p(i (i-l) = max l(i) > 

max2(i) - max {r (i~l)> 

(p,a) e B-E(i) 

If max 2(i) is zero, go to step 6. 

If max 2(i) is integer, go to step 3. 

Otherwise go to step 2. 

Step^ : Si *- max l(i) - max 2(i) + ^ 

go to step 4 * 

Step 3 i a x max l(i) - max 2(i) 

Step 4 : 'Let = (N, E(i)) 

to solve IP C G-j_5 e I] s 

(4.1) Solve HP £]G(G ^) )B 33 

(4.2) Construct an optimal solution for IP C G ( G i)? e H 

(4.3) Construct an -optimal solution U(i) for 



IP L ^i’ s i e H where 


D( i) 



Ste_g_5 ; 


Step ^j5 : 


r p q ( i ) *- max 

go to step 1. 


i r ~ 

'» pq. 


! o 
\ 



S u (3) 

3=1 P 



GHAPEB. IV 


SOLE PEOBIEMS BELA'IEP 10 P |G;f] 

4 .1 liiEROPUGIIOlT 

in tiie present chapter we identify some problems such 
as transportation problem, generalised vertex packing problem 
etc., which can also be solved using the decomposition 
techniques discussed in previous chapters. 

In section 4.2, we develop an algorithm for PjG?r~], 
when G is a tree. In this integer optimal solutions for 
decomposed subproblems of IP [Gjr j are obtained by using 
dynamic programming. An attempt is also made to calculate 
an upper bound on the number of subproblems to be solved. 

In section 4.3, a natural modification of the algorithm 

(3.4) , when the underlying graph G is bipartite, is presented. 
Por an example dual of the standard transportation problem 

is of the form IP [G,*rj on a bipartite graph G- . 

In section 4.4, we again deal with a generalized edge 
covering problem (BP), when upper bounding constraints on its 
variables have also been added to it. It is shown that HP 
can also be solved, by solving its linear relaxation and a 
weighted edge covering problem on a subgraph of G. Moreover 
linear relaxation of BP can again be solved by the algorithm 

(3.4) after making some minor modifications. 
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In section 4.5, we discuss the generalized vertex packing 
problem (GVP). Here it is shown that by a. simple linear 
transformation, this cm bo transformed to a problem of the 
form HP and thus can be solved by the algorithm modified for HP. 

In section 4.6, v/e t ake a problem whose constraints are 
mixture of constraints of generalized vertex packing and 
generalized edge covering problems . he snow that this problem 
can also be transformed to a HP . 

4.2 ALGORITHM FOR P[S$r], WHEN G IS A TREE 

Here an integer valued optimal solution of LPjGjr] is 
obtained by solving its subproblems which are generated by 
using the same decomposition technique as described in chapter 3. 
Eor solving each such subproblem, a dynamic programming type 
of algorithm is suggested which guarantees to give an integer 
valued optimal solution. Thus by solving IiP [&.$ s^e ~| finite 
number of times, an integer valued optimal solution of 
LP [G|r~| can be obtained and thus there is no need of solving a 
weighted edge covering problem. 

Since G is a tree, underlying graph of each subproblem 
will also be a loopless graph. In case, G^ is disconnected, 

LP [GMj s will get reduced to independent LP’s on these 
components and therefore can be solved seperately. Hence 
without loss of generality, we assume that G^ is connected. 
Solution U obtained by the following algorithm will 'be shown 
to be optimal and hence U(i) = U. 
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4.2.1 notations and Definitions 


Outer node 
S 1 

Ij_(p ) 

deg(p ) 

IM(p) 

f(p) 

cl 

P 


P 

4 .2 .2 
Step Q 


Step_J. 


node having degree one in G 
a subset of the node sev IT of G-p 
A maximal connected subgraph of Gn containing 
onlv the node p and nodes from the set S-, . 
Degree of p in I.., where S. is a subgraph 


of G.^ . 

{q. I (p 5 a) e g^ > 


{qjq. G s i 311(1 (P t q.) G > 

Value of the objective function for the optimal 

solution of IJP [f . (p e~] provided u = 1 , where 

e is a summation vector corresponding to the 

edge set of T^(p). 

= c + £_ c 

p P qGH(p) q 

Value of the objective function for the optimal 


solution of IiP [f ,• (p )? ej provided u. 


0 


Here c° = £_ P 

P qGH(p) 


min {c p , 


“p°> 


•?! input p i jc pJ ¥p e & ± ~J 


s 


1 


T i <- G i 


S 0 + {p jp e deg (p) = 1} 

I'ind N(p) for every node p of S Q 
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N(p) «- {qjq 6 S 1 and (p.q) e } 


P 

0 

"p 

'1 

■p 

5 1 


^ c n + 2 _ c ¥pes 
P qeN(p) 1 ■ 0 

"*■ , c n ¥ P e s 0 

q.ei'T(pj 1 
+ min (c-k c°> 


P' P 


S 1 ^ s c 

/ 


where T i -/S n is the graph T i , left after 


r o 


removing all the nodes of S Q and edges incident 
on them . 

Step_2_ (i) If ih contains exactly one node go to step 4. 

(ii) If contains exactly two nodes, go to step 3 

(iii) Otherwise go to step 1. 


Stj2jD_J3 Calculate c , c° for every p Q I. as follows 


P 
. 0 


C 4 - Z ^ C 

P qSN(p) q - 


1 


c 

qSN(p) < 1 


let the nodes of T- be p and q. 

I j-L •r- U U 

1 J-Jq 

lind V/t t = c„ + c q where lq+Lg > 1 and L 1? Lg g. {0,1}. 




p, 


o 


let W T j = min { y/ T T > 
d f2 L ± +L 2 > 1 -^2 


Set u 


P, 


1 


u 


1, 


= J { 


and go to step 5. 
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Let p Q Toe the only node of 21 . 

1 o —* 

Lind c , o and c 

P 0 p o Po 

— T 

If = c" where J e (0,1> . 

Set u = J and go to step 5. 

p ^ 

Step^^5 V/c now give values to all the other nodes of G1 . 

Select a node p which has teen given a value. 

(i) If Up = 0, set u^ = 1 ¥ node q 8 H(p) . 

(ii) If u = 1 give values to all other nodes of IT(p) 

ir 

which have not teen given values till now, as follows: 

If c^ = c^ where q e I?(p) 

Set u = 1 

q. 

Otherwise, set u^ = 0 . 

If all the nodes of have teen given values, stop \ otherwise 

again go to step 5. 

4 .2 .3 Pr oof jqf_opj^ir^ijty 

It will te sufficient to give only an outline of the proof 
for optimality of the solution obtained ty the algorithm (4.2.2)* 

Prom the description of the algorithm itself, it is easy 
to see that for any node p, 1 • ( p ) at the' time when c\ c°, c 

JP if If 

are calculated, is uniquely defined in such a way that if 
q / p, then exactly one of the following is true: 

(i) '-^U) cd T i(p) 

(ii) ' '.^(p) 

(iii) I.(q) 
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I'rom the definition 0 f c° 

P J 

that c is the optimal value of 
the stage, when the number of node 
or two . 


and o , 
p P 

^ L?i(p)J e J 

s in ^(p) 


it is obvious 
IT ow consider 
is either one 


-^ e "k contain only one node, 

Obviously I.(p ) = G. and hence c 

P n 

of IP [G . 5 e~[ . 


say it is node p o « 
is the optimal value 


CajJL 2 let 11 contain only two nodes, 
In this case, l i (p Q ) \J l.( ^ & 

Xle then give values to the nodes 


and q Q • 


P 0 and q Q , such that 

to ” 1 3114 the tot al cost for T 1 (p 0 ) (J 4 0 ) is 
in ^ c^" + c 0 ^ -i- ^ 1 


u + u 
p 

min 


*0 £V % Pc’ q 0 -p 0 


+ c 


4.2.4 fe^_^ora^a«j S iB.re A uirea to solve a subproblem 
IE[G.;e:] 


For each node q of the graph of a subproblem, we have 

to calculate c 4 , c°, 5 q . For this, only additions and 

comparisions are needed. in particular 

deg(q) additions axe necessary for computing c 1 

1 

dog(q)~l additions are necessary for computi n*r c° 

1 

One comparision is necessary for computing c 

Hence the total number of operations to compute c 1 c° c 

Q* q’ q 

for q = l,...,k, where k is cardinality of the node set of G. 

k 

<21 deg(q) ~ k+k 
q=l 

< 4(k-l) 

0(k) • 


(since G is a tree) 



G f may be dropped out because of being nonbinding in a 
subsequent problem end it may agai n reappear in a later 
problem, let El denote the number of subprobloms P r 
of A in which e^ reappears for the first time, i.c., 
o. G P r and e. ri P r ~^ for some r such that I +1 < r < 

XX* 0 

I Q + 1(10$ or for r = I +1 it appears in the problem 
I +1 

P • , Since each subproblem belonging to A, differs from 
its immediate predecessor, there is at least one edge in each, 
which reappears in it for the first time, in tho above sense , 

Hence I(k) < 2 HA 

e ^6 G- k 

V/e will first find HA for the case when, G k is 
connected. It will be easy to see from this derivation that 
the same value of can also be taken as an upperbound 

1 r 

for HA, for the case, when Gr is a forest. 

let e = (p , q ) be an edge of G k . Deletion of the edge 
0 O' 0 

e Q ■ from G-^ partitions it into two components S^(e 0 ) and 
Sg(e ). let p 0 G S 1 (e Q ) and q Q G 3 2 (e Q ). How we indicate 
some observations from the algorithm which will be used for 
the derivation of H^ . 

Here the subproblems referred to, have the subgraphs of G k 
as their underlying graphs. 

s 

(1) An edge e Q = (p 0 ,q. 0 ) gets dropped from P , 


where 
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r D g ( s 0 ) < max 1 ( s 0 +1 ) 

^ 0 1 o 


4.2 .5(a) 


where r (s ) is the requirement on the edge (p ,q ), 
1 04-0 th 

updated after s iteration. 


Inequality 4.2.5(a) is possible only if 

u p 0 (s o ) > 0 

and u (s ) > 0 

V 0 

b y 

(2) o q reappears in P 0 after being dropped out from P 

only if r (t -1) = max 1 (t ) 

1 o H o 

Since at any iteration s, (s < s < t ) 

( s— 1 ) = r ( s-2 ) - u (s-1) - u (s-1) 


P 0 1 0 


Po^o 


hence r (s-1) < max 1 ( s ) * 


Po% 

This is so because , t o ’ is the first iteration after s Q 


th 


iteration where e o reappears. Moreover if either u^ or 


u is positive in U(s-l), then 
4 0 

u (s-1) + u (s-1) > max 1 (s-1) - max 2 (s-1) 

P 0 c i 0 

Here in every iteration j, max l(j)-max 2(j) Is always an 

integer and max 1 ( i+1) = max 2(i) . Therefore edge 

s' t 

can not reappear unless somewhere between p 0 and P 0 
u (s) = u ( s ) = 0 fors ome g . 

P 0 *o 

(3) Let Gl(s ,e ) and C2(s 0 ,e o ) be two components into 
which the underlying graph of p S ° gets divided after 
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iteration# let the node p^ he in Cl(s 0 ,e o ) and lot q Q 
he in 02(s o? e o ) . 

As long as the components Cl(s 0 ,o o ) and 02(s o ,e o ) continue 

s _ 

to "bo the components of the suhpr ohlems that follow P the 
nodes p Q and q Q will get positive values because optimal 
solution for the IP on this component is obtained with the 
help of previous problem,, and hence the edge o o will not 
reappear , 

Thus if e Q reappears in t iteration, by observation 
2 at least ono edge from S^(e o ) - Cl(s ,e ) and one from 
Sg(e o ) - G2 ( s 0 j e Q ) should be added to Cl(s o} e o ) and 
C2(s o ,e o ) respectively to change the structure of these 
comp onents . 


(4) Any edge 

or q . is 
D 

problem of A. 


o. = (p-,q.) of G lc , such that either p. 

J J J J 

. • lr 

an oncer node in G“ , will be present in every 

v 

Phis is so because c.>0,¥jSG'. 

3 


low using the above properties, we will prove the 
following theorem: 


< 1 + min {jS 1 (e i )!, Is^o^l) 
where js.(e.)J = number of edges in S.(e.) for j = 1,2. 

3 1 j i 

Out l ine of t he pro of 
i-, ip 

Let P , P , . ,.,P be all the subproblems of A, 

where e. reappears for the first time. With each of the 
i o 1 

P s , s = 2,,.,,L1 we will be able to associate two edges 

f and g such that , 
s °s 
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(1) (a) f g G S^(e^) • (To) g^ g Sg(e^) 

where S^(g^), Sg(e^) ar e the components of G-^ c , obtained 

by deleting e^ from it, 

and 

(2) f g 7^ f^ D = <2, 3 j . . . , s— 1 
Sg 7^ S-j D = 2,3, . » . , s —•1 

Me will therefore be able to conclude that: 

Nj. < 1 -h min (|s 1 (e i )|, |S 2 (c.) j > , 

Prj^of: lor the sake of easiness, let us enclose some of the 

problems of A into boxes as follows: 

( i) In the first problem of each box, edge e^ reappears 
for the first time and it is also present in each problem 
contained in a box. 

(ii) In the last problem of a box, edge o. becomes non 
binding and therefore gets dropped. 

(iii) Any problem of A, which is not contained in any box, 
does not contain the edge e^ . 

let P and P denote the first and the last 

"til 

problem , respectively of the r 1 box. 

i s 

We will first show that with each P ° 9 s = 2, 3,..., IT, 
we can associate an edge f such that f 6 S-,(e.) and it 

3 S s. x 

must have been introduced in at least one problem of A, say 
in p”k, before the problem P s , and it has- played a role in 
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changing the value of node through a series of changes 

between problems P S “ 1 and P _s only. We can call such a 

change as cumulative chain effect of f on n. , It will 

s L i 

be denoted as CHE(f , t) wnere t is the iteration -where f 
was introduced and 3 , 1 < t < i., 

S — -L S 

ihis type of selection of edges will guarantee that f has 

i s 

not been assigned to any problem of A before p 3 , 

g’s can also be assigned in the same fashion. 

Note: Reappearance of an edge e^ occurs due to the change 

of the values of p. and q. both, f and g are the 

first edges which can be considered to be responsible for 

changing the structure of components containing and q_^ 

respectively, after e. was dropped out. Thus f and g 

can be considered to be responsible for the chahge in values 

irom 1 to 0 of and q^, before e^ reappears in P 

i 

Assignment of f and g to P s , (s = 2,...,N.) 

T M ... q. ^ ^ ^ w J>r .^ ^ ^ | J| ip ^ 

i 

We will first assign f to P s and then in the same 

S X 

manner g g can also be assigned to P ° . Consider any two 

boxes s^ and s^+1. Then edge becomes nonbinding in 

3 o 

O-u 

P , with respect to its optimal solution. Consider the 

3 S"K , 

components formed by the binding edges of P . Let 

01(3 , e.) and C2 M , e.) be two such components 

s b 1 s b x 

containing p^ and q^ respectively, where i> ^i^ * 

- 1 -S'h+1 

Since e^ reappears in P , it implies that between 
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s b s b +1 

p and P > at least one edge from S 1 (e.) - 01( j , e.) 

J- X 1 

has been added to 01( ,©^). let e^ = (m^,n^) be one such 

edge which has first changed the structure of Cl(j ,e.). 

s b 1 

let n^ e 0l( i ,e.), and suppose this edge e 1 was introduced 


v w *-*• ■"] V_/ ^ On 7 n / 7 ^ uu 

i s b x 

bf 

in P , where i < t^ < i 

Si -L 

0 


®b +i ■ 


Since n^ G S^(e^) 


S 1 (e 1 ) < S^e.). 


l^i- 


There are now two possible cases: 

Oase_ji; e^ is anew edge, i.e., in any problem of A, before 
t l 

P , has not been introduced. In this case take j = 1 

and go to step b. 

0_as e^ b ; 

St_ejD_ja. e^ is an old edge. By an old edge, we mean that in at 

least one problem of A, before p it has also been dropped 

out. let s., be the largest number, but less than t^, such 

S 1 

that in the optimal solution of P , e 1 ha s become nonbinding. 

Obviously s 1 < 2 , because the underlying graph- of every 

1 S-u 

subproblem between 2 and t 1 0 iterations, contain the 

s b 

component 01( 3 , e.) as a component and e^, wnich is not 

s b 1 

contained in Ol(j , e.), but is incident upon a node of 

s b 1 

Gl( 3 jcw) . let Ol(s 1 ,e 1 ) and C2(s 1 ,e 1 ) be two components 
containing m 1 and n 1 respectively and are formed from the 

S-i 

edges of P which are binding with respect to its optimal 

% 

solutions. Since e ± reappears for the first time in P , 
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it implies that between P 1 and P at least one edge 
from S^(e^) - Cl(s^,c^) which is incident upon a node of 
Ci(s 1 ,e 1 ) has been added to Ol(s 1 ,e ) . 

let e^ = (mg,ng), where ng e Cl(s^,e^) be one such edge 
which has first changed the structure of Gl(s 1 ,e 1 ) . 

Let it appear for the first time in t 111 iteration, where 

Aji 

S 1 < t 2 < t l* 

Since ng G S^(e^) 

S 2 C & 2 ) ^S 1 (ei) cl s 1 (e 1 ). 

Por eg again there may be two cases: 

Case b^ eg is a new edge , 

Case bg eg is an old edge . 

Case b^ Og is a new edge, therefore set 3 = 2 and go 

to step b. 

Case bg Let Sg be the largest number less than tg, such 
that in P e? becomes nonbinding. Now we repeat the 
procedure suggested above for an old edge e^. 

Say, after repeating this procedure 3 times, we get a 
sequence of following edges 

e l» e 2> * * * ,C 3 

Note: Por any r; 1 < r < 3 , e r = (m r ,n r ) was introduced 

t 

in P r , such that e ^ had already been dropped out, i.e., 

and n r 6 cl ( s r -l’ e r _i) * By the 0110108 of 


s r~l ^ ^r ^ ^r-l 
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these edges, it is easy to see that 

s l( e j) •*' S 1 (e i ) 4.2.5(h) 

T ,fe now claim that the above repetition of step a, will he 

done only a finite number of times. Proof is as follows: 

Since i^i( e j_)| Is 9X1 integer, from 4.2.5(b)it is easy to see 

that wo will finally get an edge e. which reappears for the 

t, 3 

first time in P 3 and for which js 1 (e.)j = <*> and 

n. e Cl(s. .,,e. -i ) . Hence m. is an outer node in . 

J J 

By above observation 4, e . reappears for the first time in 

J 

the first problem of A and will be present in every problem 

t I +1 

of A. Therefore e . is a new edge for P 3 = P ° . Thus 

J 

finally we will get a new edge e such that step a is not be 

J 

repeated any further, 
how go to step b. 


Step b W'c will first prove that the cumulative chain effect 

t . 

of the introduction of e. in P 3 can reach the node 

^ 3sp Isft+1 

only in bet ween the problems P and P . 

t„. 

In P e . has been introduced for the very first time i.e., 

t i 3 

before P no problem of A contains o... 

(i) If j = 1 

e 1 = (m^n^) where n 1 0 Cl( 3 ^ ,e . ) . 


s b x ‘ 


It is obvious that OHE (e t .) can reach the node p . 


belonging to 01(3_ ,e.) only in between P 

V b 

P D . ' Therefore, we can take f = e.» 

s b 3 


S k and 
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(ii) If 3 > 1 

t . t . 1 

(a) Consider the interval £p p , which contains 

the problems P 1 such that t. < i < t. 1 . 

3 3~x 

Sinco u is introduced in p ^ for the very first t i m a 9 

and n^ e Gl(s^ 19 e 

CHE (o^ t 3 ) -> 01(8^, o. wl ) 4.2.5(c) 

i.o., CHE (e_.,t 3 ) can effect anode of Cl(s i 1? e . 1 ) . 


H ow 


e . 
3 


6 ^l( G -j_l) o. i which is a connecting edge 


" ~ '3- t. t. , 

does not belong to any problem of £p P ^" 1 ) . 
CHE (e t 3 ) / s 2 (e ) 


Therefore 


i ,e 


3' 

CiJE (e .,t^) can not change the value of any node of 
J 


s g ( 0 j _i). 


Since node p. e S 9 (e. 1 ), 

X & J— x 

t . -j- 

CHE (e^t 3 ) / Pi in [> 3 , P ;3 ” 1 ) 

t. 1 

(b) Nov/ consider the interval [^P , 

Since n^ e Cl(s^ 2 , e j„g) 

CHE (e^, t 3 " 1 ) - Cl(s._ g , 0j _ 2 ) 

Since m.._^ is present in Gl(Sj_ 1? e j„i) 
can say, 


4.2.5(d) 



by 4 .2 ,5( c) , we 


Ci-IE(e.,t d ) -01(8^, e . m1 ) U 01(3^, e^) 

n t -4 Q 

e. o does not belong to any problem of £p , P J ~ ), and 

e. and e. 1 both are in S 1 (e. 9 ), therefore 
J D*^x x 
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Oiffi (o 3 _ 1( ti- 1 ) /s 2 (o._ g ) 

How Sg(e. - 2 )CS 2 (e.J 

Therefore CHS (e., t^) /• S 9 (e , 9 ) . 

Since is also in Sg(e^ 2 ), therefore 

OiiS (e., tp / PjL 


4.2.5(e) 


Here we can say that CHE (e t" 1 ), can not reach the node p. 


before p 


Vs 


t 9 t^ 

( c) Proceeding this way till the interval £p , P ), wo 


can conclude that 
.3 


CHE (e_. # t^) -*■ Cl(s^^j o ,i~l) ^ • * * 4/Cl(s^,e^) 


4.2.5(f) 

4.2.5(g) 


and CHE (o., t^) / p • before p 

_ ti 

(d) Finally consider the interval [_P , P ) 

Since n- e 01(3 , e .) 

r s b 1 

and m 1 8 Cl(s 1 , e^) 

By 4.2.5(f), we can say that 

CHE (e t^) - Cl( j , e ) 
j b b 

Since p is also in Cl ( 3 ,e.) 

1 s b 

CHS (e t j ) - Pi . 

Hence by 4.2.5(g) wo conclude that the cumulative chain effect 

t JL 

of the introduction of o . in P can reach the node of 
only in between P. P Sb . Hence we take 


' s-u+1 

b 


e . 
3 
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Similarly g g 's can also bo assigned to each first problem 
of a box . 

Hon co wo can say, 

l^il+min {|S 1 (o i )j, js 2 (c.)|> 

The orem 4 .2 .5 1 % 0(n 3 ) where n is the number of 

nodes in G. 

Proof: It is easy to sec that the sum 2 , min {18,(0.) I? 

. _ e er 1 

l^2^ e i^l^ w: ^- "be largest when e^’s form a graph having 

only two outer nodes* In this case, above sum will bet 

< 2 J_ 1+2+3+ ... J where k is the number of nodes 

in G-k. 


Hen cc 


i"® 2 + * * * + ^k~l - 2^+3... + (k-1) 

Therefore , 

2 (fe 2 -) (- k f 2 + 1 ) 

I(k) < — + O'- 1 ) 


= | (k-2) + (k-1) 


Since 1 « 1(1) +1(2) + ... I(n) 


n , n 

< 2 £ (k-2 ) + 2 (k-1) 

" k=l 4 k=l 

, 21 o OH 

= ~ T 2 k" + 22k - 4n1 

4= u ^_i v_i ~ J 


= 1 + 2n|n+ll ~ 4n 2 

% 0(n 3 ) 


n 



78 


4,3 i^GOttlTlll I'OR pgjf] ON A -BIPARTI 'HE GRAPrl 

Consider the problem PjGjr"], when G is a bipartite 
graph. Here a natural modif icat ion of the algorithm (3.4) 
is made, which gives an integer optimal solution for 
lp[G;r j. In tnis case subproblems are generated by tne same 
decomposition method as described in chapter 3. 

In the algorithm (3.4) at the i iteration, to get an optimal 
solution of U? jG^je], wo solve DUP £ G^) j ej , which is a 
maximum flow problem on a bipartite graph generated from G. * 
Here in this case, G^ itself is a bipartite graph, therefore 
DDPjG-.je"] is a maximum flow problem on G. . Thus the only 
modification of the algorithm is that, now wc solve DDPjjj.je"] 

_ t ‘ 1 

instead of DIP j_G(G^)jeJ , Prom this solution, as shown in 
Appendix A, a (0,1) optimal solution for IP [GA jej is generated 
from which optimal integer solution of DP [G^,s-o is found. 
Pin ally by adding these optimal solutions for subproblems, an 
integer valued optimal solution of l2?{G;r] and hence an 
optimal solution for pfc^r] is obtained. 

4 * 4 POUNDED VARIABLE JPAODDBi, ^ 

In the present section we doal with a bounded variable 
version of P [G ; r] , which is formulated as follows: 

EPnSjrl 

n 

minimize 2 c-u. 

i=l 1 1 

such that u^+ u^ > r^ ¥ (i, j) 6 G 

ui < 

u - >0, Integer 


¥ i G G 
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where , 

w . 1 s are nonnega/cive integers * 

We will first show that an optimal solution of B? \o;r] can 
also he decomposed in two parts as we ha.ve shown for an optimal 
solution of P [G;r~j in chapter 2. 


Aoove result can easily he proved hy proceeding on the 
same lines as in chapter 2 and using the following: 


Since w.'s are integers, 1/2 e will he an optimal linear 

solution for IBP [_ G.?eJ where & is a fractional component 

3 3 

of G as defined in chapter 2 for IP jGjrJj and IBP /G.^ej 
denotes the linear relaxation of HpjG^se"]. 


Alpjor it hm ( 4 . 4 ) (To solve 1 HP [& ; r~| ) 

/ v th 

As given in the algorithm (3.4), at the beginning of i . 

iteration, we will first compute gw. Let {i^,ig, , ..,i^ } he the 

set of nodes of the graph G^ of LP [Gh ; s^e ~j . We now find 

in as follows: 


+• To mr 

At the i iteration, we then solve IBP s^ej . After finding 
an optimal solution U(i) of IBP [G^;s^e] , we update r and 


min {w. , w. , s_. } 

H x 2 x p 1 


w both for all the nodes p 

p ^ 


and q of G^ * If for any node p, 


updated w becomes zero, we replace hy a large integer 

number M, so that in all future iterations u^ is forced to 
take zero value. We then give values to all the neighbouring 
nodes of p according to the updated requirements on their 
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adjoining arcs* Again after updating vi and r ¥ p,q of G, 

P P d 1 5 H J 

v;e proceed for the next iteration as before. 

Modified algorithm 4.4 gives on optimal 
solution for IBP [G^rJ . 

?,Fi?. 0 »£ • Proof is exactly on the same lines as we have given 
for lP[G;rj in chapter 3. It is easy to see that all the 
five properties P-^,Pp, . . . ,P on which wh ole proof is based, 
hold for this case also. 

4 ,5 GJIIERALIZEJJ VERIESC P ACKING. PROBLM 

Here we discuss the generalized vertex paking problem 
which can be considered to be a generalization of the vertex 
packing problem in the same way as the generalized edge covering 
problem for the edge covering problem. 

Generalized vertex packing problem is formulated as follows: 

GYP [G$r] 

n 

maximize Z c u 

p=l P p 

such that u +u ~ r pq ^ (P?l) S G 

^ > 0, Integer f p 0 G 

We now show that GYP[G;rj can be transformed to a BP[G;r , [] for 
some r' . 

let r = max {r > 

(p,q)GG P* 
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Using the transformation u = 5e-U' , Stt[S;r] transforms to 

maximize 2 c (r - u' ) 

P p' 

such that r-uj+r-u^ < r pq ¥ (p,q) e G 

U P < r Tf p e & 

Since Z c r is a constantj from the selection of r, it is 
pG& ^ 

easy to see that uf 0 ?p 6 S. G-iZP [& | can also be written 

as: 


minimize Z c u 1 
p P P 


such that 


> 2r - r pq ¥(p # q) e G 


0 < u' < r 

P 


¥ pGG 


which is same as the problem -SPjGjrj where 


PQ. 


2? ~ r pq 9 CP>^) 6 G 


Hence G-VP [&$r] can also be solved by the algorithm 4.4. 


4.6 HIXED OOIT S'JBAIK'i'S PKQHL M 


Consider the problem, 
&CW[&;rj 

minimize 


Z c u 


pGCr 


P P 


such that U p + u 4 > r pq 


u + u < r 
p q. ~ PI 


} 


°1 


u g > 0, Integer ¥ p G G- 

where and Sg are subsets of the edge set E of G- . 



v/ e snow tnat this can also be transformed to a HP p';r f ']. 
1st f = max |r pq J(p,q_) g Sgl 


Mate a linear transi ormation U = re — 1 ’ on the constraint 
set of Sg* G-CVP jG-jrJ is then equivalent to 


Minimize leu 
p P P 

such that 


11 +11 

> 

r 

p 1 


pa 

u’ +u’ 

> 

r 1 

p 1 


pa 

u' 

< 

r 

p 



u 

> 

0 

p 


u +u * 

= 

“r 


V pGG- 

V pgG such that (p,q)eS 0 for 

P P X NJr 7U/ 6 

some q * 

As in section 4 . 5 , here also we can add nonnegativity restrictions 
on u^ in P . In order to convert the equality constraints into 
inequalities, we attach a large penality cost when ever the 
equality constraint is violated. 

Thus GG'VP [g 5 rj is also equivalent to the following: 

minimize 2 c u + K f" 2 (u +u ! )-r 
pe& P P - p P P 


such that 


u -i- u > r I S., 

P 1 “ P4 j 1 

u + u > r ! / S f 

P q. - pq 

u + u 1 > r 
P P ** _ 

0 < u' < r 
“ P " 

u >0 

P 


? p e g 


which is of the form BP[&',' r 'H for some S^aph G* and a requirement 

vector r* . Therefore GGYp[&?rJ can also he solved by the algorithm 
( 4 * 4 )* 
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5 .1 BTli^nClilOT 

In this chapter we discuss computational experience 
with the algor it 1m (3,4) described in chapter three for solving 
linear relaxation of the generalised edge covering problem (P ) . 
Some observations on the behaviour of the algorithm (3.4) are 
ma.de from the computational experience gained by solving several 
randomly generated test problems. 


'Jo have shown in chapter three that in the worst possible 

case, number of linear programming subproblems generated are 

bounded by the highest requirement among all the edges of G- 

(max R = max) r . . t(i, j)eG) . Further each iteration will 

require at most 0(n ) computations where n is the number 

of nodes in the underlying graph of the original problem P . 

!hus in the worst case the computations are hounded by 

3 

l(n) £ max R. 0(n t '). In general all these generated subproblems 
are not needed to be solved separately, because they have the 
property that if for a subproblem, optimal solution of its 
immediate predecessor is feasible for it, then it will he 
optimal also. Ihus the above bound l(n) would very seldom 
be attained. Computational results have also shown that the 
average number of iterations required are usually of the 0(n) .. 
Here by number of iterations we mean number of different 
subproblems solved. 
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5 .2 OTIFOEfc PROGIL&M 

ft 

Computer program for the algorithm (3.4) is coded in 
FORTRAN -lO -DEC -lo version and is run in LSCloSO : -11 4-0 process 


system. As the main aim of the computations is to get an 
insigh b into the algorithm, the program is leapt simple and is 
not cificienb from the point of view of programming efficiency. 
As the program does not use the best data structure and other 


computer techniques, the computational time is not considered 
as an important index of performance. This program is 
structured with the help of following three main subroutine ss 


a Subroutine 'GEN' 1 

b S ubr ou t ine '1 TER 1 

c Subroutine 'VALUES’ 

a Subroutine 'GEE' It generates random test LP -problems for 

the given values of n, de and max R, where 

n = number of nodes in G 

and de = density of the graph G i .e . 

n umb er o f edge s present in G 

possTbTlT number* of~lfdges "on n’ nodes 

Random graphs G are generated with the help of pseudo 
random numbers which are distributed uniformly over [ 5 , 1 ] • 

b Subroutine J ITER' .This subroutine generates subproblems 

for all the iterations. 4 subproblem of i iteration is 
generated by computing maxl(i), max2(i) and s i as described in 



85 


chapter three for the lip problem having requirements updated 
after (i-l)^ u iteration. 

„ c , , It is used for solving a subproblem 

generated by the subroutine * I'lUR* « It first generates a 
symmetric bipartite graph corresponding to the underlying graph 
of this subproblem and then it solves the maximum flow problem 
on this bipartite graph by Ford and Fulkerson 1 s labelling 
method [i8J . 'there are other efficient algorithms available to 
solve the maximum flow problem such as, Dinic [13J , Malhotra 
[?C1 j but for the sake of program simplicity Ford and Fulkerson's 
method has been selected. 

From the optimal solution of this maximum flow problem, a 
(0, 1/2,1) valued optimal solution of the corresponding subproblcm 
is generated. Solution thus obtained is then transferred to 
the main program for updating the values of r. Us. listing 

X J 

of the program is given in Appendix B. 

5 * 3 QOI'IPU '1AH.0U ill PM WM .CE 

Performance of the algorithm is studied in terms of the 
number of subproblems solved. Information about CPU' time is 
also collected but is not considered as an important index. 
Randomly generated problems up to 900 nodes have been solved 
in reasonable amount of GPU time (3 minutes or less) . 

Main aim of the computational experimentation is to study 
the number of iterations as a function of problem parameter 
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max R, and. alts o ics variation with n and density of the 
underlying graph. 

Tesc problems are generated randomly "by varying problem 
parameters n, de and maxR. I* or each fixed value of n,de 
and max R, throe randomly generated XiP problems are solved. 
Average number of iterations ana GPU time (without problem 
generation time) are noted. 

To study the effect of max R, on the number of iterations 
required, we define the parameter : IR = maxR/m, where m is 
the number of edges in the graph, i.e., m = x de . 

I? or n = 25, 50, 100, 200, 300, random test problems having 
underlying graph of low density, such that the number, of edges 
present are n, 2n, 3n, are solved. I' or n = 50, 100, 150, 

200, problems are also solved where the graphs are of little 
higher densities, in such a way that the number of edges, present 
are £ lOn . Value of muxR for each problem (i.e. for a fixed 
value of n, m/n) is selected such that IR = .5, 1,2, 3, 4, 5 and 6 
i.e., maxR = »5m, m, 2m, 3m, 4-m, 5m and 6m. Comp u tat ion al 
results are shown in fables 5.1, 5,2, 5.3, 5.4 and 5.5. Table 
5.1 shows the number of subproblems (iterations) solved for 
various values of IR. Here it is observed that the number 
of iterations becomes almost constant when IR is increased 
beyond 2, for a specified value of n and m. 

Table 5.2 shows NR = (number of iterations/no of nodes) 
for varying values of n and max R. 
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Pig. 5.3(a), shows the variations of liR with IR. It can 
ho observed that those ratios (NR) are very close to each other 
and are almost independent of the number of nodes. Let I(n) 
denote the maximum value of NR for n nodes. This can be 
taken as number of iterations/number of nodes when 
IR = 2 f Variations of l(n) with n, for different values of m/n 
are plotted in fig. 5.3(b) It can be observed that for a fined 
m/n, I(n) is constant or is independent of n. further more 
the value of I(n) is increasing with the increase in m/n, 
but even for dense d graphs it is approximately 2, Thus we 
can conclude that the number of iterations are of 0(n) . In 
fact, such an observation can be explained by the fact that 
the number of iterations are governed by the in-differenceness 
of r. . T s and hence really are not dependent upon the magnitude 
of maxR, but on the distinct number of r.,'s in the problem. 

1 J 

To show further the behaviour of I(n) for dense graphs, 

random problems are generated for n = 50, 100, 150 and 200 

with 2 x lQ , 
n 

Table 5.3 shows the variation of NR 'with max R . for these 
problems also above observations afe confirmed. 

Table 5.4 shows average computation time in GPU( seconds) for 
randomly generated problems of size n = 25 to 800 > when m 
is taken to be 2n and max R - 2m, It is observed that even 
such large problems can be solved within a reasonable amount 
of time (less than 3 minutes) , 
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5 .4 CON PLUS I OH 

\ 

.vi'th bxio above ilieoretical and comp u. tat ional experiences 
we are now able to conclude th^ f ollowing ^regarding algorithm 
(3.4)i 

(a) I' or any II 1 , size of the first problem is very sma.ll and 

the number of nodes in successive iterations aro n on de creasing • 
luring the last few iterations subproblems on all the n-nodes 
may have to bo solved* Thor of ore on an average for the sub— 
problems 12? [CL; in gj } we can take to have n/2 nodes and 

?l§=il * U(fe(JS . 

(b) Nor different values of n and de, number of iterations 
generally varies between n and 2n for large values. Therefore 
the average number of iterations can be taken to be ^ n . 

Although the number of iterations i.e., the number of 
maximum-! low problems solved, are observed to be approximately 
, total number of generated subproblems, which depends upon 
r. . * s, may bo much larger than the number of iterations (i.e. 
the number of maximum flow problems) . Since the actual 
computational effort needed for generating a subproblem and 
comparing it with the proceeding subproblem, is too small in 
comp ar is ion to the computational effort required to solve it, 
it will not be bad to state that the complexity of the aJ-gorithm 
is bounded p olynomially by a function of the problem size 
(no of nodes) . 
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Table Ho. 5 .1 
For number of iterations 


n m/n 




__Max R/m 



_ _«J5 _ 

1 

2 

3 

4 

5 

6 

25 1 

9 

12 

15 

1 

16 

16 

17 

" 2 

13 

17 

21 

21 

22 

23 

23 

" 3 

16 

21 

23 

25 

26 

25 

26 


50 1 18 24 32 33 32 33 34 


31 37 44 44 47 46 46 


? f 

3 

35 

41 

50 

50 

49 

52 

53 

loo 

1 

37 

49 

65 

68 

68 

72 

70 

t f 

2 

54 

70 

86 

89 

86 

88 

89 

t f 

3 

68 

90 

99 

94 

104 

112 

no 

2(30 

1 

75 

105 

127 

132 

132 

135 

134 

f f 

2 

115 

145 

165 

177 

181 

178 

181 

f f 

3 

146 

169 

200 

198 

202 

215 

215 

300 

1 

116 

163 

200 

198 

209 

218 

218 

? f 

2 

171 

210 

252 

260 

264 

2 74 

272 

\ f 

3 

20 4 

256 

290 

301 

306 

305 

30 6 
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Table 5.2 

[Number of iterations/nj 


n 

m/n 

~ Y5 ~ 

'"T 

jlJax 

... 

R/rnj 

~T ' 




25 

1 

.36 

.48 

.66 

.63 

.64 

.64 

*68 

1 f 

2 

.52 

.65 

.84 

*86 

.83 

.9 2 

.92 

! t 

3 

.63 

.84 

.93 

.95 

1.02 

.99 

.93 

50 

i 

.36 

.48 

.64 

.66 

.66 

.68 

.68 

It 

2 

.53 

.77 

.89 

.83 

.94 

.92 

.92 

1 t 

3 

.69 

.81 

.99 

.99 

.96 

1.02 

1.02 

100 

1 

.3 7 

,49 

.65 

.68 

.68 

.72 

.70 

1 1 

2 

.54 

.78 

.86 

.90 

.96 

.88 

.90 

1 1 

3 

.69 

.90 

.99* 1 

.93 

1.05 

1.02 

.99 

200 

1 

.38 

.53 

.63 

.66 

.66 

.68 

.67 

t ? 

2 

.58 

.72 

.82 

.88 

.90 

.88 

.90 

T t 

3 

.72 

.84 

1.02 

.99 

1.02 

1.00 

1.03 

300 

1 

1 

! 

«i 

I 

.54 

.66 

.66 

.69 

.73 

.75 

1 1 

2 

.56 

. 70 ' 

.89 

.86 

.88 

.91 

.90 

H 

3 

.69 

.84 

.99 

.99 

1.02 

1.02 

1.02 



Table 5.3 


[lumber of itorations/nj 


n 

m/n - 


jllax R/m j 



.5 

1 

2 

3 

50 

io 

1 .1 

1,2 

1.25 

1 .4 

100 

10 

1.1 

.9 

1.29 

1.3 

150 

10 

1 .0 

1 .2 

1.24 

1 .25 

200 

10 

.9 

1.13 

1.3 

1.3 


Table 5 .4 
lor GPU time 

n_ .. J£.5? rations _ _ _ G PU time ( sec onds). 


25 

21 

.287 

50 

44 

.94 

100 

86 

3.505 

***•**»» .« •* 



200 

165 

13 .144 

300 

252 

36 .192 

400 

375 

68.52 

500 

459 

100 .9 5 

600 

49 7 

118 .23 

800 

648 

19 5.82 



Table 
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5.5 

[Number of iterations"] 


n 

density - 


Max H/m 




1 

_ 2 

3 

4 

25 

* oQ 

23 

25 

29 

30 

30 

r t 

.75 

25 

28 

36 

34 

36 

50 

• 50 

55 

60 

65 

67 

66 

i r 

.75 

54 

62 

71 

72 

70 

75 

.50 

87 

10 7 

lio 

110 

109 

! f 

.75 

87 

111 

in 

118 

‘“Tib 
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OP IIllAD COj'MUNICATIOIT Sp aLTIVING 1M1BS 


6 .1 

Hu [ 27_J has discussed tne optimal communication 
sp aiming tree problem for a given set of requirements among 
n nodes. He has also suggested an algorithm £2?£ for 
constructing such an optimal communication spanning tree. 

Here we extend the above problem to the following three cases : 

(i) Construction of an optimal communication spanning 
tree such that all tne nodes of a given set S are outer nodes. 

(ii) Construction of a new optimal communication spanning 
tree with the help of an existing one, when some of the 
requirements have been increased. 

(iii) Construction of an optimal communication spanning tree 
such that some specif ied edges are present in- the communication 
tree • 

In the last section we suggest a minor modification to 
iiu’s algorithm £27 £, for constructing a cut tree. This in 
some cases reduces, its total computational effort. 

6 * 2 liOTA'flQHS AMD DBFI NIIIOB S 

Definitions are collected from Hu's book £26 £j . 

G : Is a complete undirected graph on n nodes. These nodes 

may represent cities which need to communicate with each 
other. An edge joining nodes i and j will be denoted by 
an unordered pair (i,j). 
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"'"ij ° a nonnegative integer number assigned to every pair 

(i, j), such that r^.. = r^ where i / 3 j and i = l,2,».,,n, 
<3 ~ l 5 2 ,...,n. It represents the requirement (e.g. 
number of approximate telephone lines needed) between 
nodes i and 3. 

dj_ ^ : A nonnegative number assigned to every pair (i,j). This 

represents the distance between tne stations represented 
by the nodes i raid 3 in &, We will here assume that 
d ij = d ji = lr v ( i » 3 ) S G. 

^( r ij) ! Gra P il G witil capacity r i - on its arc (i, 3); ¥ (i, 3) e G. 

Outer : A node, whose degree is one, in the tree, 
node of 
a tree 

G(T) : Total cost of communication on the spanning tree T. 

There is a unique path from node i to node 3 in T. 

Length of this path is the sum of the distances between 

its adjacent nodes in the path. Cost of communication 

for a pair of nodes i and 3 is equal to r . . multiplied 

by the distance from i to 3 in T. Summing up these 

communication costs for all the (g) pairs, we get the 

total cost of communication on T. 

R(X,X): Is the capacity of a cut (X,X), sepe rating two nodes i^ 

and ig. Here X is a subset of N (node set of G) 

containing i^, and X is the complement of X with respect 

to IT, and R(X,X) = 2 2 _ r • , 

i 0 X 3 6 X J 

Since r. . = r.. ¥ (i,j) 6 G, 

X J J x 

R(X,X) = R(X,X). 
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f : If an edge (1,3) having capacity r . . is taken as the 

P Q. 3 

union of edges, each having capacity one in G then 
the maximum number of arcv/ise disjoint paths between the 
two nodes p and q in G will be denoted by f . 

Minimum cut sep orating p and q : 

A cut (X,X), such that p 0 X and q 0 X is called a 
minimum cut seperating nodes p and q, if R(X,X) is the 
least among the capacities of all the cuts seperating 
p and q. 

Cut : A cut tree I of a graph G(r^.) is a spanning tree 
tree J 

having nonnegative integers v- - corresponding to 

’i" J 

each of its edges (i,j), such that the followings 
are true : 

(a) : If an edge (i,j) with the corresponding number v^ 

is removed, nodes of G(r . .) get partitioned into two 
sots X and X, and 

V i3 = R(x s x). 

Also this cut (X,X) is a minimum cut seperating nodes 
i and j in 

(b) : Maximum flow f from p to q in G(r i;j ) is equal to 

min {v a , •••» v i j> * * * *’ v tq * 

where the edges of { (p ?a) i» j ) j •» •? ( t, q) } 

form the unique path from p to q in 1. 

N^(p) : (Neighbourhood of p in the tree 1). 

N T (p) = {q 0 N | 3 an edge (p,q) in 1 > . 
where N is node set of G . 
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(i) For an edge (i,j) of A, where A is a spanning 
tree on n nodes, v. . is the capacity on the arc 

-L J 

(i,j), let (X^,Ag) he a cut set separating i and 3 
in A, then 


v4 


10 


R(Il,X 2 ). 


(ii) If (i,j) $ A, v^. is the maximum flow from i 

to i in the graph A(v ) where v is taken as the 

pq pq 

capacity of the edge (p,q) in A, ¥ (p,q) & A. 


6 . 3 SOME P ' ROPBitf IE S. 01 liJE QPJPJI-.IAJj GOIgiTMI G aIIOK SPilgllJ Cj_ JDHEE3 

Hu £27"] has shown that a cut tree of G(r . . ) is an 

optimal communication spanning tree, for the set of requirements 

fr . >. Fere we will show that the converse of the above is also 
ID 

true i.e. any optimal communication spanning tree for the set 
of requirements {r^ ^ i is also a cut tree for the graph G-(r^), 
where G is a complete graph, and r.j's are taken as the 
capacities of its arcs (edges). In order to prove this, we 
first report the following observation : 


Observation (6.3.1J t Capacities of the cuts given by the 
edges of an optimal communication spanning tree are the (n-1) 
maximum flows present in the graph G(r^_^). 

Proof : let A be a cut tree for &( r j_p 911(1 let B an 

optimal communication spanning tree, let be an (n-1) 
dimensional vector whose components are the capacities of the 
cuts given by the edges of A. 
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^5 




A 

where v^ is the capacity of the cut given Toy an edge (i,j) of 
A in G(r— ) . './e will now use the following result given by 

Adolphs on and Hu £ 1 £j : 

"Given two spanning trees of a network G we shall refer 
to one as the 'red* spanning tree and the other as the 'blue* 
spanning tree, for any blue edge (p,q) there is a red edge 
(pljlf) such that (p-j_,q^) is one of the red edges which form 
the unique path from p to q in the red spanning tree. 


furthermore, a one to one mapping f^ among the blue 
edges and rod edges can be established which satisfies the 
above condition". 


Using the above mapping f f 1 from the set of edges of A 
to the set of edges of 33, we can construct an (n-1) dimensional 
vector Vp whore a component of Yg is the capacity of the cut 
given by the edge f 1 (i,j) of 33. Ye claim that the vectors Y 1 
and Vg are the same. We prove it as follows : Since A and 33 
both arc optimal communication spanning trees, total costs of 
communication on A and B are equal. Hu £27 £ has shown that 
the total cost of communication on a spanning tree is equal to 
the sum of the cut-capacities given by its (n-1) edges .Therefore 

e7 i = eV 2 
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where e is the (n— 1) —dimensional summation vector* Suppose, 
a component; of Yg is greater than the corresponding component 
of V l5 then there must exist a component of Yg which is strictly 
smaller than its corresponding component in Y 1 . let this 


component of be y"; i.e. 

r l r 2 


B 

Vf ( r l» r 2 



V 


A 


r l r 2 


6.3(a) 


where f(r 1 ,rg) = P , <1 1 such that (p ! ,q» ) e f 1 (r lJ r 2 ) . Hdge 

fl( r l , rp) Xi the Xrom r^ to r 2 in 3, therefore it gives 

a cut set (X,X) sepe rating nodes r^ and r g in G(r^) and 

B 

capacity of this cut is vt, v. Since A is a cut tree for 

i^r^,ig; 

G(r. •), therefore vf’ „ is the capacity of the minimum cut 
10 r i r 2 

sepe rating nodes r^ and r 2 in G(r i ^). Hence the inequality 
6.3(a) gives a contradiction. Thus Y^ and Yg are not different. 


Since A is a cut tree, components of V 1 are all the (n-1) 
maximum flows present in the graph G(r. •). Henoe the edges of 

X J 

any optimal communication spanning tree give all the (n-1) 
maximum f 1 ows for G ( r i ^ ) . 

theorem XArAA). : Any optimal communication spanning tree B 

for the set of requirements ^r^, is a cut tree for the associated 
graph G(r i; p . 

Proof : By the observation (6.3.1) edges of .13 correspond to 

the (n-1) minimum cuts seperating all the (g) pair of nodes of 

G(r. .)• Hence an edge of B represents a minimum cut seperating 
x 3 

a certain pair of nodes of G(r^) . 
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I'low Toy the theorem (9.2) [26 3 proved by Hu, this spanning 
tree 3 will have all the properties of a cut tree. 


Hence 3 is a cut tree for G(r. - 

v ID 


). 


Li 


6.4 OP f II i A1 COMhUMCA'iTO! SPiL'lIITO- THESE HAVILG All THE HOIElS 

DF X : spa oiPinfT seX" W its‘ otW Yarns 


we now present an algorithm for constructing a spanning 
tree such that it contains all the nodes of a specified subset 
S of N, as its outer nodes, and the total cost of communica.tion 
on it, is also minimum among such spanning trees. 

Let S = { s 1? Sg, . . . ,s lc ,s k+1 , . . .jS-^ } be this specified set. 

P or a spanning tree B(p) of G(r . .), where p is a integer number. 

x 3 


S Q (p) = { s I S q G s and it is an outer node of B(p)>. 
P or a node s ^of S-S 0 (p), define 

NB B(p) (s n .)= {Cl e o |(s ,,a) e B(p) and q. $ S (p) > . 
P+1 P+1 

Let B(0 ) = A, where A is a cut tree for (^r^), 
and let, S Q (0) = (s k+1 j * * *> s b wilere |S J = b 


Starting with an optimal communicat ion spanning tree B(O) for 
{r . . > we will construct a sequence of spanning trees 

13 

B(l), . . .,B(k) such that 


s 0 (p) C s 0 ( p+1 )’ 
|S 0 (p+l)| = |s 0 (p) I + 1 
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and j S 0 ( lv ) j = b and S 0 (k) = S, 


For eirery p, p - 1, ...,k , 'spanning tree B(p) will be obtained 

from tiie spanning tree B(p~l) by the following algorithm : 

Algor it hm _( 6_.4) s 

: led G S and it is not an outer node of B(p-l). 


Find = max 

^ je NB^P- 1 ) 

Step.l : Let j e i B (p _1 )( s ) 

the edge (s^,j) from B(p-l) . 


( , v 8 ^- 1 ) l 


(s p K 


) -p 




6.4(a) 


Add an edge (j,q ) and delete 


After each repetition of step one, we get a spanning tree 
and after s ) | repetitions, a spanning tree, in which 

the node s^ is also an outer node, is obtained. We will denote 
this spanning tree by B(p). 

By definition 1JB B ^'~ 1 ^( s n ) f\ S 0 (p~l) = and 

A.' «> 

a 0 s ), it is easy to see that, 

P P 


s 0 (p-i)e:s 0 (p) 

and |S 0 (p) | = |s 0 (p-l) | + 1. 


For example let B(p-l) be, 





Figure 6 .4(1) 



lol 


where I-^Ig and Ig denote the parts of B(p-l) . 
Let 

S 0 (p-1) = {SfjSg, . . * s s p4 ,s k+1 , . .. } s b > 
NB B (P~l)( Sp ) = { q pj i l? 3l> 


and v 


B(p-l) 

S P 4 P 


max > v 


U^P- 1 ) ,B(p~l) _B(p— 1 


s„ 1-, 5 s.. 


lYfe ' Vl ~p 

Then by the above algorithm (6.4), B(p) will be 


IT ow 



S Q (p) = { s-^Sg j • • •) s p„i» s p j s ]£+l» • * * > s b 


Pj?0Qerty~JJ)s4j : Prom the construction of B(p), it is easy to 
see that for every edge (i,j) of B(p) which was also in B(p-l) 
Dut not incident upon s^ there, following is true : 


r B(p) = 


6.4(b) 


3-D 3-D 

eorem (6.4,1) : Por any two nodes p and q. such that neither 

or q belongs to S (}c), value of the minimum cut given by the 
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tree B(k) is also a .minimum out; so per at ing nodes p and q 
in G(r^) . ■ 


Proof : We pro-ve it by induction. Let 3(o) = A vine re A is 

an optimal communication spanning tree, for G(r..). 3 y 

1 1 ) 

theorem (6.3.1) A will be a cut tree for {r. Therefore 

1 3 

above theorem is true for B(i), when i = 0. Let the theorem 
be true for B(i-l), i > 1 we will now show that it will be 
true for B(i) also. 


Select arbitrary two nodes of G, say p and q, such 
that neither p nor q is in 8 (i). 

Let (p , i^, ig, . . ,,i ,q) be the unique path from p to q in B(i-l) 
and let s^ be the node which is made an outer node in B(i). 
Consider the following cases : 


Case_ 1 : s ^ does not lie on this path. 

: s i d 063 on 'fcbis path. 

Case, 1 : By construction of B(i) from B(i-l), it is easy to 
see that the unique path from p to q in B(i) is the same, as 
it was in B( i-1) . 



By 6.4(b) it can also be written as 

= min S 

L p 1 
= ^ i " 1 h 


. . . , V. 


B(i-l) 


i r q. 
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Since neither d nor q iq in q ( i 'i A ,, 

y 4. xs in S 0 (i-1) and the Theorem is assumed 

to he proved for k > i~l v B(i-l) .• + , „ Li 

* - ’ v pq I s xne capacity of the 

minimum cut seperating p and q. in G(r. .). Hence v B(i ) is also 

p q 

the capacity of the minimum cut seperating nodes p and q. in 
G ( r i^‘ 

: fet q_^ he the node adjacent to s i in B(i). Here 

again two more cases are possible : 

( i) does not lie on this path. 

(ii) q. i does lie on this path. 

JLCA), • dn B( i~l), let the path from p to q he 

(p, i-j_, . . . , i Q , s^, j Q , . . ., j r ,q) . 


For example let 



In B(i) path from p to q will be (p,^, . . .,i 0 ,q i , j Q , . . j r ,q) . 
Therefore 



6 


# 
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B CO 



By construction of B( i) from B(i-l), it is also easy to see that, 


Y 


\ ( i) = 3( i-l) 


Vi 


Vi 


v B(i) = B( i-l) 
^o^i j 0 s i * 


Hence hy 6.4(h) we can say, 

y B(i) = v B(i«l) # 

pq pq 

Since the theorem is assumed to he true for B(i-l) and 

neither p nor q is in SHi~l), and therefore is 

r x o v/7 pq_ pq 

the capacity of the minimum cut seperating nodes p and q in 

G ( r i j) * 

Cas e (i j) : In B(i-l), let the path from p to q he 

(p,i 1} ...,i o ,s i ,q i ,j 0 , ...,3 r »4)* Path from p to q in B(i) will 

therefore he (p , q, . . i Q , q i , i Q , * . . , j r ,q) • 

For example if B(i-l) : 
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Then £(i) will he 



i3y the choice of q. (from equation 6.4(a)) we know that 


°i^i 1 o s i 

Hence can also he written as, 

n n •* 


PI 


B( i-1) 
PI 


mm 


(J(i-l) B( i-l) B(i-l) J( 1-1)1 

vpq Vi ’Vo ’ v J 
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i3y construction of B(i) from B(l-l) it i s easy to see that 

v B(i) = B(i-l) 

1 o q i v i 0 s i ’ 

By 6 ,4 (Td) ? cut: capacities for the other common edges of the 
two paths from p to q in B(i— 1) and B(i) respectively, are 
the same . hence 


B(i-l) = p3(i) 

p q. pi 

Hence v^^ is the capacity of the minimum cut seperating nodes 
p and q in G(r . 

Since the above theorem has 'been proved for an arbitrary 
i, Therefore it will be true for i.= k also. 


□ 


Wo will now show that B(k) is the required optimal 
communication spanning tree. 


: Total cost of communication on the spanning 

tree B(k) is less than or equal to the total cost of communication 
over any other spanning tree of G(r . .) which has all the nodes 

X J 

of S as its outer nodes. 


?£P.°,£ ! Proof is by contradiction. Let C Q be a spanning 
tree having all the nodes of S as outer nodes and let the total 
cost of communication on 0 Q be strictly less than the total 
cost of communication on £(k). Por every e S, let t.. and v^ 

be its neighbours in C 0 and B(k) respectively. 

how consider the trees B* and 0*, where 



and 


10? 


J3> = B(k) 


U 


s i e s 




°h 


°o - 


VJ {(s^tpi 


S i e s 


i.e. B‘ is obtained from B(k) by removing all the edges (s^,v^) 
where s- G S. Similarly 01 is obtained from 0 o by removing 

X 0 0 

from it all the edges (s-,t.) where s- 6 S. 


We now use Adolphs on’s ^ 1 Z] mapping 1 f^' from the set 
of edges of B* to set of edges of C’ . II ode sets of B and 
are the same. Bor every edge (p.,p.) of B' , there is an edge 
f-,(p.,p.) = ( q. ,q.) of C’, which is in the unique path from p. 
to p^ in C^. f^_(Pj_jPj) also lies' in the unique path from p^ 
to p^ in C Q . Therefore the edge f-j_(p^ } p^) = (q^jq^) in CM 
gives a cut seperating p. and p. in G(r . .). By theorem (6.4.1) 

v n^n^ ^s ca P ac ity of the minimum cut seperating p - and p . 

P iP j 13 

in G(r ±;j ). 

c 

Therefore v^) < Vg° ¥ (p^p.) G B* 

0 

Hence 2 v®^) < 2 v 0 

(Pi»Pj) e 3' P i P 3 ~ (q i? qp e G' Q ^i^ 


How for any node s^ of S 5 


V 


B(k) = 

s i T i p G G °.i J 


2 r 


s,P 


s * t . 

^ *1 n 


i i p G G- °i J 


S-iP 


and 
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Therefore 


Thus 


Z = 2 


0 . 


Y_ 


S,V, _ „ n 


s i e s °i*i S i 6 s 


x x 


(s i ,t i )G C 0 


S v B ( fc ) + 2 < 2 + 

(p i5 p.) 6 3' p i P D s i S S s i v i " (s i ,s.)e G' 0 


c. 


, s iS . 


+ 2 v ° + 

S-j^ 0 S S i i 


and hence 


t b « < 


V, 


0 

(Pi»Pj) e B(lc) ' p fPj " (s^sp 0 C 0 ' S i S 3 


6 *4( c ) 


Hu £2 7^] has shown that the total cost of communication over 
a spanning tree is sum of the capacities of cuts given by its 
edges. Thus the ine qualitjr 6.4(c) gives a contradiction. Hence 
B(k) is the required optimal communication spanning tree. 


□ 


Ei x amp le_ (6.4) * This example is borrowed from Hu £27[J 



Big. 6 ,4( 7) 


(Edges hawing zero capacities are not shown). 

Optimal communication spanning tree A obtained by Hu for the set 
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of requirements {r- . > is 

r 3 



Let the specified set S whose nodes are required to be outer 
nodes in the communication spanning tree, be = {2,5,6 >. 

S o (0) = <f> . 

Let us first select node 2 to malce it an outer node i.e. 




= 18 


Hence = 1. 

Therefore optimal communication spanning tree having node 2 
as an outer node, is 



B( 1) 

Eig . 6 .4(9 ) 





S 0 (l) = (2 > 

node 5 of S is not in S (1) 
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Therefore q.g = 3 . 

Hence following is the optimal communication spanning tree 
having node 2 and node 3 as' its outer nodes : 



Hence B(3) is 





S 0 (3) = (2,5,6) = S. 

Thus B(3) is the required optimal communication spanning tree. 

6.5 OPTIMAL CXDMMUNICATIOIT SPAITOi'G THEE POR TjHB UPRATED SET 

OF TE quiWeMs ' * 

In this se ct i on an alg or i t iim f or c on s tr u ct ing an op t imal 
communication spanning tree for the updated set of requirements 
{r \ .} with the help of an optimal communication spanning tree 

X J 

A for the old set of requirements {r. is presented. 

Suppose the requirements {r. .} have been updated to 

o a q 

as follows : 

(1) Requirement between nodes p^ and pg is increased by 5^ 
(6q is a positive integer). 

(2) Requirement between nodes q 1 and q g is increased by 

. (6 g is a positive -integer) . 

(r) Requirement between nodes z 1 and z p is increased by 6 

-A* JO 

(6 r is a positive integer). 


Let B 0 = {(p 1 ,p 2 )> (iq? lg), • • •} ( z 1? Zg) > be the set of 
all the pair of nodes for which requirements have been increased. 
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Me will first find out those parts of the tree A, which can 
he made to he present in an optimal communication spanning tree 
for the updated set of requirements {r^>. 

let N(B Q ) = {i | (i,j) S E 0 for some 3 > 

Thus N(E 0 ) = {Pjl, p 2 > q.i» ci 2 , . . .jZi,z 2 >. 

let T he the smallest connected subgraph of the given 

p+1 

optimal communication spanning tree A which contains all the 
nodes of N(E ) . let be the set of all the edges of A which 
are incident upon a node of I hut are not present in T^ + ^. 
Deletion of all the edges of B 1 partitions the spanning tree A, 
into several components say T^,Tg, . . .,T + i# let (t^,p^) he an 
edge of A which joins T^ and I ^ where t^ 6 and p^ 6 ^p + i» 

Set can also he written as ; 


B 


“1 = Ug.pp !(% P p 6 A where tj. 6 p t 6 I p+1 > . 


It is easy to sec that if a minimum cut (X^,Xg) sep era-ting 
nodes i^ and jp in G(r.j.) does not contain any edge of E 0 , then 
the cut (X^,Xg) will he a minimum cut seperating nodes i^ and 
in G(r! .) also. Using this we will now show that an optimal 
communication spanning tree for {r! . > can he constructed, which 

»L J 

contains all the edges of Tp \] Tg\J ... T . lor constructing 
a cut tree for G(r! . ) as suggested by Hu £ 26 ? we 0321 start 

with an arbitrary pair of nodes, let us first take the pair 
ti,Pi« Edge (t n -,p i ) of A gives a minimum cut ((T. >, {1^}) 




seperating nodes t- and p • in G(r..), where { 1 ■ > denotes the set 

=X X X J X 
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of all the nodes contained in T^ and { T^ > is the set of all 
the nodes of G which are not in [ 1^ . This cut does not contain 
any edge of B Q , therefore it will be a minimum cut seperating t^ 
and p • in G(r.- .) also. After repeating it for every pair of 
nodes t^ and p ^ of B^, structure of the cut tree constructed 
so far, for G(r ! .)_ will look like : 



/ r n \ 
( 1 / 



T ) 
VP / 


l ig * 6 .5(1) 


Gut given by any edge (i 1? i 0 ) of \J T. does not contain any 

r ^ i=i 1 

edge of B Q , therefore the same cut will be a minimum cut 

seperating nodes jL and i Q in G(r’- .) also. Hence all the edges 

& x 3 

,p 

of U T. can be put in a cut tree for G(rl .). 


low to obtain the complete cut tree for G(r'. •), we have 

3-3 


to add edges in | r -^p + ij on ly * Inis will be done by finding minimum 
cuts seperating different pair of nodes contained in {Tp + ^>. 
These cuts will be obtained in a simpler graph G Q of G, which 
is obtained by condensing all the nodes of a (i / P+1) to a 
big node. Using the above results, a stepwise algorithm for 
constructing the required optimal communication spanning tree 
is given below : 
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Ste p Q : Find i* 1 the given optimal communication spanning 

tree A. Find the minimum cut (X^,Xg) seperating any two nodes 

contained in I l1f in the graoh G . 

p+i 7 , ° - o 

S tep 1 : Select a component, say X-^ containing more than one 
node of I Find a minimum cut (Y^,Yg) separating two nodes 

of CD which are contained in X^» Low do branching from X^ 
or Xg according as XgC_ or Xgd_ Yg £26 £j . For example : 
if Xg(£ Y^, structure of the cut tree will look like 





Fig. 6.5(2) 


Repeat step 1, till each big node of the present cut tree 
contains not more than one node of CO^^. Since all the cuts 
are determined in a simpler graph G Q , all the nodes of a 
will be present in one big node of the present cut tree, 

Step^jE : Let a big node X Q of the present cut tree, contain 

a node p of 1 - and all the nodes of { T . I i J {1. > \ j ... 
r P +J - 3j_ u 3g 'W 

>. Since we have already shown that ({CD.}, {CD.}) is a 
D-jj 3 3 

minimum cut seperating t- and p - in G(r'. .) also, therefore t. 

j J J D 

of T. will be connected to this node p of 1 where t. is 
j i r p *r ± 3 
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a node of I. such that the edge (t 


3^3 


) is in A for some 


fq of T p+1- 

Thus an optimal communication spanning tree can be 
obtained by finding minimum cuts for (r-1) pair of nodes only, 
where r is the cardinality of the node set of T 


PjDroljLarj_£6^._52 : If all the edges of E q are the edges of A, 
then A itself will be an optimal communication spanning tree 
for the updated set of requirements {r? •}» 


6 .6 


OPTIMAL 
BLGBSf ’ ~ 


COMHUli I CATION SPAELTILCf TREE 1IAYIHG SOME P5ESPBCIPIB.D 


In the present section, we deal with the problem of 
constructing an optimal communication spanning tree, when it is 
required that, some specified pair of nodes should be adjacent. 
We will further assume that these specified pair of nodes do 
not form any cycle. Construction of such a spanning tree may 
be needed due to one of the following reasons % 

(1) Requirement between such specified pair of nodes are of 
higher priority, 

(2) Informations to ho communicated be tween such pair of nodes 
are classified and therefore they should not pass through any 
other node . 

let B^ = { ( i 9 3 ) | i and j are required to be adjacent > . 

We will show that a out tree ’I 1 obtained for the graph 
G(r ! .) is the required optimal communication spanning tree, 

X J 

where 
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r 


id 


r. + ( Z r..) x (n~l) , if (i,j) 6 E- 

3 (i,j) S G- 13 1 



otherwise . 


First we will show that this cut tree contains all the edges 
of E^ . 

Theorem (6.6,1) : Every edge (i,j) of is present in T. 

Proof s let A be an arbitrary spanning tree on n nodes and 
having all the edges of let B be a spanning tree on n 

nodes, not having atleast one edge of E^. let such an edge be 
(i ,j 0 ). length of the path from i to in B is atleast two. 
Therefore the cost 'of communication from i Q to j on B is 
atleast 

= 2 x ( Z r.,) x (n~l))+ £r. i . 
i? 3 3 ° J o 

Total cost of communication over B is greater than 


( Pi I + 1) x .( Z r ± , X (n- 1)) + 

i 5 3 

where R^ is the sum of the costs of communication on B with 
respect to the requirements {r^}. 

Total cost of communication on A with respect to the set of 
requirements { r . } is 

d 

l®ll ( Z r, x (n-l))+Hg 

i,3 3 

where Rg is the sum of the costs of communication over A with 
respect to the set of requirements {r. . > , It is easy to see 

X J 

that 



Since > o, total cost of communication over B is greater 

than the total cost of conmunication over A, with respect to 
the set of requirements { r ! . > . Since 1 is a spanning 
tree having minimum cost of communication for the set of 
requirements {r 1 . . > , it must have all the edges of E. 

X j 

life will now show that this cut tree 1 for G(r^) 
is the required optimal communication spanning tree. 

HI 

‘.Theorem 6 . 6 .2 . A cut tree 1 for the graph G(r^) is o;f 

least total cost of communication, with respect to the set of 

requirements {r. . > , among the spanning trees of G, and 

j-3 

containing all the edges of E^ . 

Projsf . Take an optimal communication spanning tree A for 
{r. . which contains all the edges of E 1 * Now again we 
use Adolphs on’s mapping f f-,' j~ l[] from the set of edges 
of 1 to the set of edges of A. 

Obviously f 1 ( i, 3 ) = ( i, j) T (i,j) 6 2^ 

Every edge (i^ig) of - j such that (i^ig) 4 E i gives 
a cut separating nodes i^ and ig in G(r^) also. This 
cut may or may not be a minimum cut seperating i^ and ig 
in G(r. . ). Erom the definition of {r*. . >, it is easy to 
see that the capacity of this cut in G(r. .) is less than 
or equal to the capacity of any other cut in G(in^), seperating 
i^ and ig and not containing any edge of E^. Cut given by 



118 


the edge f^(i, 3 ) of A, for every edge ( 1, j ) of I, such 
that ( i, j) ^ % does not contain any edge of S^. Let us 
denote the capacity of a cut given by an edge (p,q) of CO in 
G ( r ii ) "by V Hence 

^PH = V PH ¥ 6 2 SUCh that ( p ’ q ) $ E l 


And hence 


such that 


~ T < 2 

(p j q) 8 f P ^ 

(p,q) 4 E i 


v. 


f(p 


',4') 


where $ 


f (p'»q.‘) = p’q’i (p'jq 1 ) = f 1 (p>q.)« 


It is easy to see that in G(r! .) cut given by any 

X J 

edge (i^,j^) of 1 where ( i^, j 1 ) e will be a minimum 
cut seperating i 1 and X in G(r. .) also, and hence 

J- J- X- J 


- 1 r 1 „ 

v. . = f v. . — 2 r . . 

X 1 3 1 ‘ 1 1 3 1 ( i, j )SG 13 

value of the minimum cut seperating 


(n-l) 3 

i^ and 


will be the 


Therefore 2 

(PjCl)ef 

(p 1 ,q' ) = f 1 (p,q), 


” I 1 

51 ' (p,r,)eA ^(P'.4') »her 0 

which gives a contradiction to the 


optimality of A. 


Hence T is the required optimal communication spanning 
tree . -j 

6.7 constru cting 

a^urF ’tree] 

Hu £26]] has given an algorithm for constructing a cut 
tree of a graph G(r..)« In’ this algorithm (n-l) maximum flow 

^ J 
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problems are needed to be solved, where n is the number 
of nodes in the graph Of, In this section we present some 
results which in some cases can be used for reducing the 
number of maximum flow problems required to be solved for 
constructing a cut tree . 

We also specify an order in which nodes are to be selected 
for solving the corresponding maximum flow problems. Ihis way 
of selection in some cases reduces the number of maximum flow 
problems to be solved. Since the selection of nodes in this 
specific order does not require any extra computational effort, 
Hu's ["26 3 algorithm can bo modified by this proposed order. 

[Theorem (6.7.1) Any node i 1 of H, such that, 


2 r. , = max ) 1 r . , 

jeG i 3 leer [jee 13 


Will be an inner node (i.u. its degree is greater than ono) 
in at least one optimal communication spanning tree for {r^}. 

proof. Let Tq be an optimal communication spanning tree 
for the sot of requirements \ r . . V such that iq is its outer 

V. 

node. Let (iq,i ? ) 6 Iq 


o 


'( -*-2 ■ 
V, S' 



fig . 6 .7( 1) 


from !q, let us construct a spanning tree lg, just by 
interchanging the positions of iq and ig • 
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Ai U Ag = {i 1 ,ig >= { d I 36G and neither 3 = i^ nor 3 = ig > 


Total cost of communication over T^ is 

°( T l) = r i li2 + r i„A,+ X i,A, + r i,A, + r i 0 A. + r - ' + 


■ i 1 A 1 • 1-jAg 


’H A i 


-2^1 i l A 2 

r W • 


Total cost of communication over Tg is 


C('lp) = i + r. A + r i A + r i a + r i A + r i A 
^ 1 1 2 X 2^1 X 2 A 2 1 2 a 1 X 2 a 2 At 


+ 


r iiA g + 'h 


where for any node i of G 


r . 


iA< 


= 2 r . . 


jGA 1 


ID 


and 


iA r 


DGA c 


r . . 
ID 


where, Z^£Zg[] is the sum of the cost of communication over 
all those edges of T^|^Tg^ which are incident neither upon 
i^ nor ig. from the construction of Tg it is easy to see 
that Z^ = Zg * Now there may he two cases: 


(i) 


( ii) 


Z r . • = 

3 GG 43 

z 

DGG 

r V 

2 r. - > 
DGG H* 

z 

DGG 

r ig 5 
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pasAijLJ, 

Since 

2 r . . 

3GG H 2 


and 

2 r - 
jSG 2 J 

• 

» • 

0('A) 

= c(T 2 ) 


-12 HH X 1 A 2 


0-0 ” -L • i i _L - 

1 2 Vl Vl 


Thus Tg is an optimal communication spanning tree whore i^ 
is an inner node . 

Gane__ J'jLiJ^ Since T^ is an optimal communication spanning 
tree 0(1^) < O(Tg) . It implies 

r . • + r . , + r . . + r . . + r . , + r . , + r . . + 0 1 

H 2 11 1 A 2 1 2 ii l 1 2 a 2 HH Hz 1 

< r. • + r. . + r. + r. + r. . + r. ( . + r. , + 0 o 

IrjAg ^" 1^1 3-p) XgiLc) 

r . - + r . . + r . ,, < r . . + r . ,, + r. . 

. . H X Z HH HH ” HH HH HH 


2 r . • < Sr..* 

36G H* " DSG ^ 


This contradicts the fact that 2 r . . > I ri . . 

3GG X 1 J 30 G- x 2 3 

Hence node i^ will be an inner node in at least one optimal 

communication spanning tree for fr > . 

13 

AAAAXLAA.AAL ^ or 321 y optimal communication spanning tree 
with an edge ( i 0 * D 0 ) such that node i is an outer node; 
following is true 


2 r, , < 2 r . v . 

k£G 1 o A k£G- 

Pr_oof . Follows immediately from the theorem. 
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The above result can be used for constructing an optimal 
communication spanning tree as follows : 

If in the process of construction of a cut tree, we got 
an outer ’big -node* containing only two nodes say i and 3 , 
then there is no need of finding a minimum cut sen orating i 
and 3 in G-(r..), This is so because, now by comparing 

X J 

2 r.-, and 2 r we can decide which will bo the outer 

keGr ^ keG 

node with the other as its neighbour in the final cut tree. 

In iiu's algorithm [2 7], no preference is given for the 
order in which a pair of nodes should be selected for determining 
minimum cuts. It is quite possible that we may first get a 
minimum cut ( { i o > , { 1 ” }) i.o., first we get an outer node 

i and then after solving one more maximum flow problem, wo 
got Its neighbouring node with the help of the minimum cut 
( ( i Q5 3 0 > > { } )' 

(Thus two maximum flow problems had to be solved to determine an 
outer node and the corresponding edge. 

It would certainly be better if some how wo first get 
( { i 0 j 3 0 > » { r o VJ 0 '>) instead of ({i Q >, {T’ > ) because with 
this cut, now by using the theorem (6.7.1), outer node can be 
determined without solving any more maximum flow problem. 

Thus we can say that the number of maximum flow problems 
to be solved would quite often be reduced if ’would be’ outer 
nodes are not branched out before their adjacent nodes. For 
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this we suggest the following selection procedure for branching 
from a big-node, or in the begining itself, To introduce 
edges between the nodes of a set N we select two nodes 


o ? °o 


from it as follows; 


1 r . . 

36& x o 3 

2 r . . 

jeo V 


max 

feiL 


= max 

ieli 


o 


l D'e 

j 2 

\jec- 


2 r . 
a ID 


5 ] 


r . . 2 r . . ^ 2 r. 

13 36 G 13 jGG X o 3 


Next we will show that in some cases, the theorem (6,7.1) 
together with the above selection procedure of nodes can 
reduce the number of maximum flow problems to be solved. 


Suppose in the process of the construction of a cut tree 
we get a cut (A^,A^), where = {i-^igji^ > , and R^ =R(A 1 ,A 1 ) 


and 


1^ = {%>) 



2 

3G& 




2 

dgg 



2 

je& 




2 

DSA^ 



such that R. R. <r_ R- * 

X 1 ~~ x 2 x 3 

With the help of the theorem (6.7.1) it can be shown that 

R. > Rr . This can be proved as follows: One node of A^ 
1 3 ” A 1 

say i 3 is going to be adjacent to the big node A 1 (i.e., 
tree formed by the nodes of a 1 will be incident upon i 3 ) . 
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Now changing the positions of i v and. A- a better communication 

O ± 

sp anning tree can be obtained, which is a contradiction to the 



(c) 



R i 


1 


C ase_ . Since i^ is the only node of 1^ for which 

R. > R , therefore this will be the only node of A^ which 
1 3 ~ 

can be adjacent to the big node . 


Again with the help of the theorem (6.7.1) we can say that the 
cut tree will take either of the following structures * 

W * 


(i) 


(ii) 







fig. 6.7(4) 
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Since for the "big node A = ( 9 i 99 i„} , It. < R. < R. , by the 

-L & tj 1q lr^ 

above suggested selection procedure , minimum cut seperating i v 
and ig will be determined first. As we know that this 
minimum cut is going to be present in the final cut tree if we 
arbitrary select nodes from {i^i^f,} for finding maximum flows, 
therefore by the above structures of the cut trees (i), (ii), we 
can say that the minimum cut seperating ig and ig in 
G(r..) will either be ({in}, {A^, or ( {i^, ig},{ A^, i^}) . 

from both the above cubs, big node can be opened without 

solving any other maximum flow problem, 'thus in case (a) a 
big-node containing three nodes is opened up by solving one 
maximum flow problem instead of two. 


Case b Since R. > A 

x 2 " 'A 1 

and R. > R 

x 3 ~ A x 


Either of the nodes ig and ig can be adjacent to A^. Row 
again using the theorem (6.7.1) we can say that following are 
the only possible forms that the nodes of a. can take in the 




(ii) 



J-A/U 



Pig. 6,7(5) 

Therefore the minimum cut seperating nodes i„ and i 9 

O (C 

(selected according to the above selection procedure) in 


G(r . .) 
v iD y 

will either 

be 

(i) 

( f A.^, 

ig 5 -“-g ^ , {ij_ 1 ) 

( ii) 

( A^j 

V ? { 1 3 5 1 1 > ) 

(iii) 

( { , 

ig } ? { j ig ^ ) 


In all the above three cuts, big node A^ can be opened up 
without solving any more maximum flow problem. Thus in this 
case also, big node A^ is branched out by solving only one 
maximum flow problem instead of two. 


Case ( c) In this ca,se, since 
can be adjacent to the big node 


R. > R , 
X 1 ~ \ 


an-? node of a* 


A^ . Using the theorem (6,7.1) 


we can say that the nodes of A^ can take any of the following 




( iii) 

( iv) 



Tiv * 6.7 ( f,\ 





Therefore the minimum cut sepera,ting nodes ig 
either he 


12 7 


and ig will 


(i) 

( { A^ p ig } $ 

{ il? 

i 3 * ^ 

(ii) 

( ^ * ig ^ 5 

{ i X 5 

i 2 > ) 

or (iii) 

( { ig ^ ? { A^ $ 

1* 5 if-r 

> ) 


In the cut sets (i) and (ii), node A^ can he opened up 
without any extra computation. But for the cut set (iii), 
decomposition of can he either of the form ( i) or ( iv) . 

Thus in this case only, we have to check, whether 

li( { i-^ i f { }) ^ E. ( { ig , ig } , { A^ , i-^ } ) 

If in the hegining , nodes of G are renumbered as 

i 1 , ir, ... i. such that R. > R. > ... R. } the above 
s. 6 n 1 i ” x 2 “ 11 

selection procedure does not need any extra computational 
eff ort . 

Hence in general, Hu’s algorithm £26 j can he modified by 
selecting nodes in the above specified manner. 

Examp,!©.. 6^7 



fig * 6 .?( 7) 
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6 

R. = 2 r . . 

1 d=l 13 

R 1 = 18 j Rg = 19, R 3 = 15, R 4 = 14 , R g = 18, R g = 20. 

Arranging the nodes such that their R. 1 s become in the non 
increasing order 

6,2,1, 5,3,4 

vie therefore first find maximum flow from node 6 to node 2 in 
G(r..). Minimum cut seperating 6 and 2 is ( {1,2 }, {6, 5,4,3 }) 

J- J 

having capacity 17. 'therefore the starting structure of the 
cut tree is 



Since R ? > R^ node 2 can not be an outer node in the cut tree, 
therefore we can open the big node {1,2 } 

Nov/ we solve a maximum flow problem selecting nodes 6 and 5. 
Minimum cut seperating these nodes in G(r. .) is ( {6,2,1 > , 

**“ ti 

{5,3,4 >) which has the capacity 15. 

Therefore form of the present cut tree is 




Rig . 6 .7(10 ) 
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or 

where 
S ince 




node 4- can neither he adjacent to A, nor it can be adjacent 
to node 3. Therefore the only possibility* - to open (5,4 1- is. 




Thus with the above selection procedure of nodes, we could 
get an optimal communication spanning tree (cut tree) only 
by solving three maximum flow problems, instead of 5. 
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AEPEMIi A 


mEC G ( s i)i e 1 which is a maximum-flow problem on a 
bipartite graph, can be solved by any of the efficient algorithm 
of Karp £ 3o_J» fulkerson £18^] and Malhotra £33 . fere we 

show that an optimal solution of its dual i.e. of UP L ! e J 

can be constructed directly from an optimal solution X of 
DU? £ G( je even if X is not an extreme point solution, 
following is the procedure to generate an optimal solution of 
LP [[G(G^);e 3] from an optimal solution of DU £&(U) je 23 : 


Algo r it hm ( A . 1 ) ; Let B'^U Kg denote the node set of the 
bipartite graph G(G^). 


0 


Step 1 : 


Set u^ = 0 


1£ 2 X P1 < °D ¥ P 6 U 1 


qSlb 


p 


Set Up = 0 if 2 x < c ¥ q e Kg 

! pe»i 

(i) If Up = 0 and there is an edge (p,q) 6 G(G^) such 

o 

that q has not been assigned a value, set u^ = 1. 

2 - 

If u^ = 0 and there is on edge (p,q) 6 G(G^) 

1 

such that p has not been assigned a value, set u = 1, 

1 2 ^ 

(ii) If u^ = 1 and > 0, set u^ = 0, provided node 

q has not been assigned any value. 

2 1 
If u^ = 1 and x > 0, set u^ = o provided 

node p has not been assigned any value. 


Repeat step 1, till no more nodes can be 
assigned values. If all the nodes have been assigned 



1 2 

values, stop | -vector (U ,U ) with the components 
1 2 

Up, u^ for p e H-l and q. e ilp gives a (0,1) valued 

optimal solution for 13? £ G(G^) ? *e j| . Otherwise go to 

step 2. 

1 

Step 2 : Set u = 0 for every node of I'j 1 , which has not been 

ass igne d any value . 

9 

Set u = 1 ? Q e ii 0 ) if there is an edge 
q. 2 ’ 

(p,q) S G(G^) such that p g Nq and has been assigned 
zero value, but q G Kg has not been. 

2 

Set u = 0 for all the nodes of fil to which 
0. <c 

values have not been assigned till now. 


1 2 

Obviously all the components of (U ,U ) are non-- 

1 2 

negative. By the construction of the vector (U ,U ), there are 

1 2 

'only two possibilities ° } for (U ,TJ ) : 

( i) It is infeasible for IP LG(Gq) s e 13 % 


Infeasibility can occur only if some or more constraints 

are violated, i.e. there exists at least one edge (p 0 ,q. 0 ) 

— 12 
in G(G-) such that u =0 and u = 0. 

^ o ^-o 

(ii) It is not optimal for UP £G(Gq) ;e : 

By construction of the vector (U^,U^) it is ensured that 
if the primal has a positive surplus 


(i.e.c - 2_ x > 0) for a constraint, then the 

q P G G-CG-p) 

corresponding dual variable is zero. Hence in case 
(U 1 ,!! 2 ) is not optimal but is feasible for Iip£G(G-^); e [] , 



there will exist at least one constraint in the dual 


with the positive surplus such that tne corresponding 

primal variable is also positive, let (p ,q 0 ) be this 

1 2 

constraint such that u = 1, u. =1 and x >0. 

P 0 % Po^o 

We will now show that at any step of the algorithm (A.l) 
feasibility condition (i) and complementary slackness condition 
(ii) are not violated. 

Step^o i Since X represents a maximum flow on the bipartite 
graph Gr(Gr^), we can not have any edge such that both the nodes 
on which it is incident, are unsaturated. Hence there is no 
edge (p ,q ) such that 




(c) 


<1 


L 

e 




G(S a ) 


P 0 4 


c and 2 „ x 
p o p G G(G i ) pq o 





In case (a), flow X can be improved by sending 


min { c 

*• v* 


L „ 


X 




q e G(G i ) p o cl q o p e G(G X ) ~ pq 


x nn } along the 


edge (p 0 ,q o ) which contradicts the optimality of X. 


In case (b) one of the nodes from p' and q Q would have 
been assigned the zero value in step 0 and therefore the other 
node will be assigned value one in step 1. Hence again none of 
the constraints can be violated. 


In case (c), since 2 _ x n - c and 

q e G(G i ) p o ci p o 


2 „ 


x„ 


p e G(G.) pq o q 


c , neither p nor q will be assigned value 


zero in step o« From step 1, since u~ = 0 ? there must exist 

q o 


an edge (p x ,q ) in G(G.) such that x > o and u - " =1. How 

P 1 q o 

1 9 

= 1 implies that there is an edge (p 1 ,q 1 ) such that u = o. 
p i r r q . 


1 

i 

Pi 


u 


In case 2 „ x v < c we stop 
k G G(G. ) q l 


If 2 _ x. _ = c we continue as before till we obtain a 

k e GCcy icq i q i 

node q for which 2 x < c . 

S p 6 G(G.) pq s q s 

Due to the finiteness of the graph it is certain that we 

will be able to get such an unsaturated node q . 

s 



If I - x = c . there must exist an edge (p-uq.,) 
q G G(G.) p l q p l 11 

2 

having positive flow and u =1. So w we proceed as before. 

q l 

Thus due to the finiteness of the graph G(G^), we will be able 

to get a node p such that i _ x _ < c 
' S q. G &(G i ) p s q p s 

Similarly starting from node p , v/e can get a node 

such that Z x _t < c r. Thus a path of odd length from 

p G G(G i ) pq t q t 

p g to q^ can be found such that both the nodes p g and q_^_ are 
unsaturated w.r.t. X. Thus as before, flow X can be improved by 
alternatively increasing and decreasing flows by a suitable 
G > 0, along the edges of this path. This will contradict the 
optimality of X. 

Step 2 : By the description of the algorithm itself, v/e can 

say that both the nodes p Q /and q Q of an edge (p 0 ,q Q ) G G(G^) can 
not get zero values in step 2. Also for any edge (p 0 ,q Q ) G G(G^) 
such that p Q or q Q gets value zero either in step o or 1, other 
node must also get a value one in step 1 only. 

Again from the algorithm, it is easy to see that for any 

edge (p 0 ,q 0 ) such that x > 0 5 either both the nodes have 

p o q o 

been assigned values till the end of step one, or both remained 
unassigned. 

Step 2 assigns values zero to all the nodes of which 
have not been assigned any value. Thus case (ii) will not occur 
in this step also. 



Since the algorithm (A.l) is such that after assigning a 

value to a node in step 'i 1 it will not change its value in the 

. 12 . 

following steps, therefore we can say that (TJ ,U ) thus obtained 
is feasible for U? $e 2] and also satisfies all the 

complementary slackness conditions with respect to X. 

Hence (U 1 ,!] 2 ) is an optimal solution for UP []G(G- i )?e ] . 
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